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Preface 


The  aim  of  this  study  was  the  measurement  of  neutrons  from  the 
high-energy  x-ray  machine  used  for  cancer  therapy  at  the  Wright- 
Patterson  Air  Force  Base  Medical  Center.  Thin  gold  foils  inserted  into 
flux  integrators  (polyethylene  cylinders)  were  used  to  measure  the 
neutron  output  from  the  x-ray  machine  and  the  results  reported.  To 
ensure  the  validity  of  these  neutron  measurements,  a  series  of  calibra¬ 
tion  experiments  were  done  on  various  parts  of  the  detection  and  count¬ 
ing  system;  these  are  also  reported. 

This  document  is  the  culmination  of  a  long  struggle  on  the  part  of 
myself,  the  author*  Of  course  I  wish  tc  acknowledge  the  contributions 
by  my  major  professor.  Dr.  George  John.  The  idea  for  this  study  was 
his,  and  I  am  grateful  for  his  forbearance  in  dealing  with  this  recalci¬ 
trant  student.  I  wish  also  to  thank  Capt  Ken  Wohlt,  medical  physicist 
at  the  WPAFB  Medical  Center,  without  whose  help  I  could  never  have  done 
the  actual  neutron  measurements  on  the  x-ray  machine.  Let  me  not  forget 
also  the  friends  who  gave  me  encouragement  to  complete  this  work*  Most 
importantly  I  thank  all  those  family,  friends,  peers,  teachers,  preach¬ 
ers,  counselors,  and  others  who  believed  in  me  during  my  formative  years 
and  beyond.  They  gave  me  the  confidence  and  strength  to  believe  in 
myself,  so  that  I  was  not  totally  discouraged  by  the  skepticism  of  oth¬ 
ers.  I  dedicate  this  effort  to  those  who  saw  the  potential  in  me  and 
encouraged  me  to  live  up  to  it.  I  also  dedicate  it  to  all  who  struggle 
through  adversity,  yet  persevere  to  finish  the  task*  To  them  belongs 
the  final  victory  as  only  they  can  know  it. 

Richard  A.  Jackson 
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Abstract 


The  goale  of  this  study  were  (1)  to  verify  calibration  factors 
determined  by  previous  experiment  for  relating  counts  from  activated 
gold  foils  to  neutron  fluences,  and  (2)  to  use  the  method  to  determine 
neutron  fluences  under  various  operating  conditions  of  the  Varian  Clinac 
1800  (18-MeV)  electron  accelerator  at  the  WPAFB  Medical  Center.  For 
comparison  similar  measurements  were  made  of  the  CGM  Saturne  1  at  Ket¬ 
tering  Memorial  Hospital,  Dayton,  Ohio. 

Cadmium-difference  pairs  in  a  standard  graphite  pile  were  employee 
to  verify  the  counting  efficiency  of  the  Geiger-Mueller  beta  counting 
system  used  to  count  the  foil  activity.  They  were  also  used  to  deter¬ 
mine  the  relation  between  the  thickness  of  the  gold  foils  and  the  activ¬ 
ity  induced  in  the  foil  (self  shielding  and  self  absorption). 

A  PuBe  neutron  source  was  used  to  calibrate  the  response  of  flux 
integrators  in  which  gold  foils  are  activated  by  exposure  to  neutrons. 
The  response  of  these  integrators,  which  are  cadmium-covered  polyethyl¬ 
ene  cylinders,  15  cm  in  height  and  diameter,  is  not  strongly  dependent 
on  the  neutron  energy  between  20  keV  and  20  MeV.  The*  experimentally 
determined  calibration  factor  is  2970  ±  180  n-mg/cm^  for  converting 
foil  saturation  activity  per  mg  to  neutron  flux. 

The  calibrated  flux  integrators  were  used  to  measure  neutron 
fluences  from  the  accelerators  at  both  medical  facilities.  With  a  20  x 
20  cm  field  size  and  100  cm  target  distance,  the  Varian  Clinac  1600 
fluences  ranged  from  ( 1.80±0« ll)xl0^  n/cm^-radx  (4.29i0.27  mrem^/rad^) 
at  20  cm  from  the  beam  axis  to  (8.96±0. 56)xl04  n/cm^-radx  (2.13±0.14 
mremn/radx)  at  100  cm  from  the  beam  axis.  These  were  measured  on  the 
treatment  table  and  without  patients.  Beside  the  treatment  table  (about 
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35  cm  from  the  beam  axis)  the  detectors  averaged  ( 1* 235±0«010)xl0^ 
n/cm^-radx  (2.938±0*024  mremn/rad^)-  Experiments  at  the  table  side 
location  during  patient  therapy  were  of  the  same  order  of  magnitude  as 
those  done  without  patients,  lending  credibility  to  the  technique. 

For  the  10  MV  Saturne  1,  an  average  fluence  of  (5.985±0.035)xl0^ 
n/cm2-radx  (2.497±0.014  mremj^/radx)  was  found  at  20  cm  axial  distance 
on  the  table,  20  x  20  cm  field  size.  At  100  cm  axial  distance  with  the 
same  field  size,  an  average  fluence  of  (2.720±0.048)xl0^  n/cm^-radx 
(1.134±0.020  mremfi/radx)  was  found.  These  values  are  somewhat  greater 
than  half  the  fluence  from  the  Varian  Clinac  1800  for  the  same  condi¬ 
tions* 
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NEUTRON  FLUENCE  AND  DOSE  FROM 


A  VARIAN  CLINAC  1800  ACCELERATOR 


I  Introduction 


ProblM 

The  problem  addressed  in  this  thesis  study  is  the  determination  of 
the  neutron  flux  and/or  dose  output  from  an  electron  accelerator.  Spe¬ 
cifically  the  goal  was  to  measure  the  neutron  dose  a  cancer  patient 
receives  while  undergoing  x-ray  treatment  with  the  high-energy  electron 
accelerator  at  the  WPAFB  Medical  Center. 

Activation  of  gold  foils  was  the  method  used  for  these  neutron  mea¬ 
surements.  When  gold  foils  are  exposed  to  a  neutron  flux,  such  as  from 
the  electron  accelerator,  some  of  the  gold  atoms  absorb  a  neutron  and 
become  radioactive.  Thus  it  is  necessary  to  relate  the  foil's  activity 
(as  measured  by  some  specified  detector  system)  to  the  neutron  fluence 
that  induced  said  activity. 

The  solution  of  the  problem  above  was  easily  divided  into  four  main 
steps.  First  a  series  of  gold  foils  were  activated  in  a  neutron  pile 
with  a  known  thermal  neutron  flux.  This  activity  was  used  to  determine 
the  efficiency  of  the  detector  used  to  count  the  foil  activity.  Second 
was  a  calibration  of  the  response  of  the  neutron  flux  moderators  to  a 
known  neutron  flux  (not  the  thermal  flux  in  part  one  above).  Third,  the 
calibrated  neutron  flux  moderator  was  used  to  measure  the  unknown  neu¬ 
tron  flux  from  the  Medical  Center's  accelerator.  Last,  a  neutron 
f luence-to-neutron-dose  conversion  was  applied  to  give  the  expected  dose 
a  patient  would  receive.  These  steps  are  explained  in  detail  later. 


1 


Background 

The  WPAFB  Medical  Center  has  a  high-energy  electron  accelerator  - 
specifically,  a  Varian  Clinac  1800  -  which  is  used  in  radiation  therapy 
for  cancer  patients.  To  produce  x  rays  the  electron  beam  is  targeted 
onto  a  material  which  has  a  high  atomic  number,  typically  copper  or 
tungsten.  The  Varian  Clinac  1800  has  a  tungsten  target.  Rapid  deceler¬ 
ation  of  the  electrons  in  this  target  produces  bremsstrahlung  x-ray  pho¬ 
tons.  These  high-energy  photons  irradiate  the  cancerous  tissue  as  part 
of  the  therapy. 

To  direct  the  x  rays  into  a  beam,  a  filter-collimator  assembly  is 
used  to  channel  the  x  rays.  These  assemblies  are  usually  of  lead  or 
tungsten.  The  Varian  Clinac  1800  uses  tungsten. 

If  the  energy  of  the  photons  exceeds  the  binding  energy  of  the  neu¬ 
trons  in  the  target  material  then  neutrons  can  be  expelled  due  to  photo- 
nuclear  reactions.  Most  nuclei  have  binding  energies  between  7  and  11 
MeV  (NCRP  31:14).  Lead  and  tungsten  -  both  used  in  accelerator  heads  - 
have  binding  energies  near  7  MeV  (McGinley  et  al,  1976:397);  the  lead 
cross-section  peaks  at  13  MeV  (NCRP  31,  1964:14).  Some  neutrons  are 
also  produced  by  electron  interaction,  but  the  probability  for  ejection 
of  neutrons  directly  by  this  method  is  about  100  times  smaller  than  for 
photons  (NCRP  31,  1964:14).  The  WPAFB  accelerator  produces  electron 
energies  of  either  6  or  18  MeV  for  x-ray  treatment.  Since  most  cancer 
therapy  is  conducted  at  18  MeV,  neutrons  are  produced.  This  is  an 
undesirable,  though  unavoidable,  side  effect.  These  neutrons  cannot  be 
collimated  as  the  x  rays  can;  thus  the  patient  receives  a  whole-body 
dose  of  neutrons.  It  is  desl cable  to  know  the  dose  of  neutrons  received 
and  to  make  certain  that  it  remains  a  negligible  part  of  the  total  dose. 

To  measure  the  neutron  dose,  knowledge  of  both  the  neutron  flux  and 


2 


the  energy  spectrum  of  the  flux  are  needed.  A  detection  system  for 
measuring  neutron  flux  in  an  accelerator  environment  must  meet  two  con¬ 
ditions.  First  it  must  be  insensitive  to  the  intense  photon  flux  pres¬ 
ent.  Secondly  it  must  be  sensitive  to  the  neutrons.  Neutron  detection 
by  activation  of  gold  and  indium  foils  meets  both  conditions.  Gold  and 
indium  have  a  low  cross  section  for  photonucleai  reactions  and  a  high 
cross  section  for  neutron  absorption;  thus  both  materials  will  absorb 
neutrons  while  being  unaffected  by  photons.  The  maximum  activity  that 
can  be  induced  in  the  foil  is  equal  to  the  rate  of  radionuclide  produc¬ 
tion  by  neutron  absorption;  this  maximum  activity  is  called  the  satura¬ 
tion  activity.  The  rate  of  neutron  absorption  is  directly  proportional 
to  the  flux  of  neutrons  across  the  foil;  the  activity  is  proportional  to 
the  fluence  (integrated  flux)  of  neutrons  over  a  given  time.  Thus  the 
induced  foil  activity  can  be  used  to  calculate  the  neutron  fluence  when 
certain  conditions  prevail.  If  the  spectrum  of  neutrons  is  known,  the 
whole-body  neutron  dose  can  be  determined.  Combining  this  fluence  with 
estimates  of  neutron  interaction  with  tissue  allows  determination  of  a 
whole  body  neutron  dose. 

A  prior  researcher  at  AFIT  working  on  the  problem  of  neutron  dose 
from  medical  accelerators  used  Bonner  spheres  (polyethylene  spheres  of 
various  diameters)  to  determine  the  neutron  spectrum,  and  flux  integra¬ 
tors  (cadmium  covered  polyethylene  cylinders,  15  cm  in  height  and  diame¬ 
ter)  and  cadmium-difference  pairs  to  measure  neutron  flux  (Rossano, 
1989:23).  Though  unable  to  measure  the  neutron  fluence  at  the  base 
medical  center,  he  did  measurements  at  Miami  Valley  Hospital  to  test  the 
neutron  measurement  methods  developed.  He  reported  dose  rates  of  1.5 
millirem  per  rad  of  photons  (Rossano,  1989:55). 
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ScQpa  of  Thesis 


Part  of  the  goal  of  this  thesis  is  to  verify  and  improve  upon,  if 
possible,  the  factors  found  earlier  by  Rossano  to  relate  counts  from 
activated  gold  foils  to  neutron  fluences  (which  caused  the  activation). 
This  requires  using  much  of  the  same  equipment  and  procedures  that  he 
used  for  his  report.  The  other  goal  is  to  use  the  flux  integrators  to 
measure  the  neutron  fluence  at  the  WPAFB  cancer^therapy  accelerator. 

To  deduce  the  flux  of  neutrons  that  activated  a  particular  foil  the 
neutron  spectrum  must  be  known  or  certain  assumptions  must  be  made  about 
it.  In  this  continued  work  the  spectrum  was  not  to  be  determined. 
Instead  the  assumption  was  made  that  the  photoneutron  spectrum  is  simi¬ 
lar  to  that  of  fission  neutrons  (McGinley  et  al,  1976:397;  Nath  et  al, 
1984:233;  Rossano,  1989:18-19).  The  average  energy  of  such  photoneutron 
spectrums  is  around  1.0  MeV  (Nath  et  al,  1984:233;  Holt  et  al, 

1979:428).  Also  the  detector  system  used  has  a  uniform  response  over 
the  neutron  spectrum  present.  Specifically  the  flux  integrators  "flat¬ 
ten  out"  the  spectrum  and  so  lessen  the  dependence  of  the  foil  activa¬ 
tion  on  the  spectrum  present  (McGinley  et  al,  1976:399;  Bruninx, 
1970:658).  The  spectrum  determined  previously  by  Rossano  was  assumed 
typical  of  what  will  be  encountered,  with  an  average  energy  of  approxi¬ 
mately  1  MeV  for  the  neutrons  from  the  electron  accelerator  (Rossano, 
1989:52) . 

The  flux  integrators,  gold  foils  and  cadmium  covers  used  in  this 
project  are  the  same  ones  used  by  Rossano.  This  equipment  is  described 
in  later  parts  of  this  paper,  but  their  usage  is  briefly  described  here. 

The  bare-  and  cadmium-covered  foils  were  exposed  to  neutrons  in  the 
standard  graphite  pile  in  the  basement  of  Building  470.  Since  the  ther¬ 
mal  flux  in  the  pile  is  )cnown,  one  can  calculate  the  saturation  activity 
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of  these  pairs.  The  known  activity  can  then  be  used  to  determine  the 
counting  efficiency  of  the  Geiger-Mueller  detector  system  used  for 
counting  activated  gold  foils. 

The  flux  integrators  containing  gold  foils  were  exposed  to  a  known 
neutron  flux  from  a  bare  PuBe  source  in  the  dome  area  of  Building  470. 
The  activity  of  these  foils,  determined  by  the  calibrated  Geiger-Mueller 
counter,  was  used  to  obtain  the  calibration  constant  for  the  flux  inte¬ 
grators. 

Measurements  at  the  hospital  were  made  with  and  without  patients 
being  present.  These  measurements  were  used  to  calculate  the  fluence  of 
neutrons  under  various  operating  conditions  of  the  accelerator,  from 
which  a  dose  can  be  calculated  if  the  spectrum  of  neutrons  is  known. 

Outline  Thesis  Pocuaent 

The  remainder  of  this  report  gives  the  theory  and  methods  used  for 
confirming  and  applying  foil  activation  for  neutron  monitoring  of  the 
electron  accelerator.  Chapter  II  summarizes  the  theories  behind  foil 
activation  for  determining  neutron  flux,  and  introduces  conversion  fac¬ 
tors  for  calculating  neutron  doses.  Chapter  III  describes  the  equipment 
and  experimental  methods  used.  Chapter  IV  gives  the  results  of  the 
neutron  measurements.  Chapter  V  sums  up  the  conclusions  from  the 
results,  along  with  suggestions  for  further  investigation.  A  series  of 
appendices  follow. 

Uncertainties  in  data  values  were  determined  by  propagation  of 
error  (Bevington,  1964:56-64)  except  where  otherwise  noted.  A  different 
method  was  used  only  in  one  circumstance,  where  the  application  of  prop- 
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agation  of  error  produced  unreasonable  results  in  the  uncertainties.  In 


that  instance 
uncertainties 


a  statistical  formula  was  used  (Bevington,  1964:93). 
are  given  at  the  one  sigma  level. 


All 
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II  Theory 


Neutron  Activation  Method 

For  this  report,  neutron  fluxes  and  fluences  are  measured  by  foil 
activation.  This  technique  is  useful  for  a  variety  of  situations  in 
which  other  methods  are  impractical.  For  instance,  the  environment  to 
be  measured  may  be  hazardous  to  humans  for  the  time  of  measurement,  or 
it  may  contain  other  types  of  radiation  that  would  flood  a  real-time 
detector.  An  electron  accelerator  falls  in  the  latter  category. 

Gold  has  excellent  properties  for  detecting  thermal  neutrons  by 
activation.  It  has  a  single  stable  isotope,  ^^'^Au,  that  has  a  large 
absorption  cross  section  of  97  barns  for  thermal  neutrons  and  a  low 
cross  section  for  photonuclear  reactions.  The  activation  product, 

^^®Au,  decays  with  a  half  life  of  2.6935  days.  In  99.975%  of  its 
decays,  ^^^Au  decays  by  beta  emission  to  an  excited  state  of  mercury 
(NCRP  58,  1985:489).  In  about  95%  of  the  decays  the  excited  mercury 
atom  emits  a  411-keV  gamma  ray;  in  about  5%  an  internal  conversion  elec¬ 
tron  of  about  340  keV  is  produced.  These  modes  of  decay  make  it  rela¬ 
tively  simple  to  measure  activation  by  measuring  either  beta  or  gamma 
decay. 

Whatever  material  is  used  in  neutron  activation,  the  theory  is  the 
same.  The  reaction  rate  R  for  neutron  interactions  in  a  volume  V  is 
given  by  the  neutron  flux  4)  times  the  total  activation  cross  section 

or 


(1) 


(Knoll,  1979:766).  The  activation  cross  section  in  this  equation  is  the 
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average  cross  section  for  the  entire  neutron  spectrum  since  a  depends  on 
the  neutron  energy.  One  can  also  express  this  as  Z-A^o  where  N  is 
the  total  number  of  target  atoms  per  volume  in  the  material  being  acti¬ 
vated  and  a  is  the  flux-averaged  microscopic  cross  section.  Given  the 
cross  section  and  fluxes  are  energy  dependent  terms  then  we  must  express 
the  reaction  rate  as  an  integral  over  energy, 

R-Nj  f  o(f)(J)(f)df  (2) 

Jo 

(Price,  1964:314)  where  Nj  represents  the  total  number  of  atoms  in  the 
target  material. 

The  microscopic  neutron  cross  section  for  gold  has  a  1/v  (where  v 
is  velocity)  dependence  for  thermal  neutrons,  that  is,  neutrons  less 
than  about  0.5  eV.  Thus  av  K  where  K  is  a  constant*  It  is  advanta¬ 
geous  then  to  choose  some  reference  cross  section  and  velocity  cTo.i^o  r 
and  express  the  cross  section  in  terms  of  these.  The  0  subscript 
indicates  the  most  probable  value  for  a  Maxwellian  distribution.  When 
the  integral  for  reaction  rate  is  divided  into  a  thermal  and  epithermal 
region^  the  reaction  rate  term  in  the  thermal  region  becomes 

— <l>(£)d£  (3) 

Jo  V 

r^-^ooVo 

=  Nr  -^niE)vclE 

Jo  V 

ro.5 

=  NtOoVo  J  n{E)dE 

Jo 

^  NjaoVon,k 
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where  nth#  number  of  neutrons  per  unit  volume  with  energies  less 

than  or  ecjual  to  0.5  eV,  is  called  the  thermal  neutron  density. 

Now  the  material  being  irradiated  with  neutrons  also  undergoes 
radioactive  decay  as  the  radioactive  atoms  are  formed.  The  time  rate  of 
change  of  radioactive  atoms  is  simply  the  rate  of  activation  minus  the 
rate  of  decay,  or 


dN 

dt 


R-KN 


(4) 


where  N  is  the  number  of  radioactive  atoms  and  K  is  their  decay  con¬ 
stant.  Assuming  N=0  at  time  t»0  the  solution  to  this  equation  is 


yV(0  =  ^(l-e'"') 


(5) 


The  activity  A(t)  is  then 

/l(0  =  ^A/(0  =  /?(l (6) 

After  a  very  long  time  -  say,  10  half-lives  or  more  -  we  have 
for  the  saturation  activity.  This  value  will  be  used  later  in  calculat¬ 
ing  the  neutron  flux.  It  is  impractical  to  irradiate  the  material  to 
its  saturation  activity,  but  this  value  can  easily  be  calculated.  If 
the  time  of  exposure  ,  the  time  between  end  of  irradiation  and  begin¬ 
ning  of  counting  ,  and  the  time  of  counting  in  the  detector  tc  along 
with  the  net  counts  (total  counts  minus  background)  are  known  then  the 
saturation  activity  may  be  calculated  according  to  the  relation. 
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A  cadmium  cover  of  about  0«02  inches  thickness  serves  to  block  thermal 
neutrons  while  allowing  epithermal  and  fast  neutrons  to  irradiate  gold 
foils  within  the  cover.  By  irradiating  gold  foils  with  and  without 
cadmium  covers  one  can  measure  the  difference  in  activity  induced  by 
only  epithermal  neutrons  and  activity  induced  by  all  neutrons.  This 
difference  represents  the  thermal  neutron  contribution.  In  practice  a 
correction  factor  must  be  applied  to  the  epithermally  induced  activ¬ 
ity  since  the  cadmium  does  absorb  some  of  the  epithermal  neutrons.  This 
correction  factor  is  dependent  upon  the  thickness  of  the  cadmium  cover, 
and  to  a  lesser  extent  on  the  thickness  of  the  gold  foil  (Price, 
1964:337).  Values  reported  in  the  literature  range  from  1.013  to  1.04 
for  cadmium  covers  20  to  40  mil  thick  over  gold  foils  up  to  about  4  mil 
thick  (Price,  1964:337;  Martin,  1955:53;  WAOD,  1960:21).  Previous 
researchers  have  used  Fcd*^*04  when  working  with  the  AFIT  neutron  pile 
(Rossano,  1989:63;  WAOD,  1960:21)  and  cadmium  covers  of  nominally  30  mil 
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(Rossano,  1989:63).  Thus  by  replacing  and  rearranging  in  the  above 
expression  the  induced  thermal  neutron  activity  may  be  written  according 
to 


A'j-Oot/on,h  »  /l.-Fcd/l.(Cd) 


(9) 


where  ^.(Cd)  is  the  saturation  activity  induced  in  the  cadmium  covered 
foil.  At  this  point  a  correction  Fth  Also  be  made  because  of  the 
neutron  flux  depression  due  to  absorption  by  the  foil.  However,  the 
contribution  by  thermal  neutrons  to  the  activity  of  20-mil  cadmium- 
covered  gold  foils  is  less  than  0.1%  of  the  total  activity  (WADD, 
1960:21).  This  uncertainty  is  negligible  in  comparison  to  other 
uncertainties  in  this  experiment,  so  it  may  be  neglected. 

From  Eq.  (9)  the  thermal  neutron  density  may  be  written  as 


n 


ih 


FcdA^jCd) 

N  jO qVq 


(10) 


Now  the  total  thermal  flux  is  given  by  .  Since  the  average 

velocity  in  a  Maxwell-Boltzmann  distribution  is  1.128  times  the  mc-t 
probable  velocity  (i.e.  1.128t;o  )r  the  previous  equation  may  be 

rearranged  to  give  (with  F^h  included  for  completeness) 


1.128F,, 


A^-FcqA.jCd) 

NtOq 


(H) 


Using  activity  per  unit  mass  is  a  bit  more  convenient  than  entering 
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the  total  number  of  atoms  in  each  foil.  The  total  number  may  be 
expressed  in  terms  of  mass  by 


Nr 


MassNa 


(12) 


where  is  Avogadro’s  number,  AW  is  the  atomic  weight  of  gold  and  Mass 
is  the  weight  of  the  gold  foil.  This  substitution  into  (11)  gives 


IA28AW 


A^-Fc^A^jCd) 

OoN  aMciSS 


(13) 


By  introducing  a  specific  activity,  5.  such  that 


5. 


A. 

Mass 


(14) 


Equation  (13)  may  be  rewritten  as 


4*  (A 


1 .128/11/ 


S.-fcd5,(Cd) 

OoN  a 


(15) 


This  is  the  eq[uation  used  in  this  study  for  the  standard  pile  foils. 
The  known  for  the  pile  is  used  to  calculate  5.  for  use  in  the  cali¬ 
bration  of  the  Geiger-Mueller  counter. 
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Mcutron  Modcratorg  and  ThcrMal  Mautron  Meaguraientg 


Neutrons  emitted  by  a  neutron  source  are  often  divided  into  two  or 
three  groups,  depending  upon  their  energy*  Less  energetic  neutrons  are 
usually  designated  as  thermal  neutrons.  More  energetic  neutrons  are 
designated  as  epithermal  or  as  fast*  The  dividing  point  between  thermal 
and  epithermal  is  generally  taken  at  about  0.5  eV;  at  this  energy  the 
absorption  cross  section  for  cadmium,  which  is  often  used  to  block  ther¬ 
mal  neutrons,  has  a  very  steep  drop. 

Gold  is  useful  as  a  thermal  neutron  detector  because  of  its  high 
cross  section  for  thermal  neutron  capture  (approximately  97  barns).  In 
a  fast  neutron  environment,  the  neutrons  must  be  slowed  down  if  the 
detector  depends  on  thermal  neutron  capture  cross  sections.  To  slow 
down  (or  moderate)  the  neutrons,  a  material  is  needed  whose  atoms  are 
approximately  the  same  mass  as  the  neutron.  If  the  scattering  particles 
are  of  a  different  mass  than  the  neutron,  near-inelastic  scatters  occur 
whereby  the  neutrons  retain  their  high  energies.  Most  organic  materials 
are  adequate  moderators  since  hydrogen  makes  up  a  large  portion  of  cheir 
atomic  composition.  Polyethylene  is  often  used  as  a  neutron  moderator; 
indeed,  the  flux  integrators  used  for  this  project  have  a  cadmium  cover 
over  a  polyethylene  interior.  The  cadmium  blocks  thermal  neutrons  so 
essentially  only  fast  neutrons  enter.  Once  in  the  polyethylene,  the 
neutrons  collide  with  hydrogen  atoms.  Energy  exchange  to  the  hydrogen 
slows  the  fast  neutrons  to  thermal  energies.  The  resulting  thermal 
neutrons  are  isotropic  in  direction,  due  to  scatter  after  giving  up 
energy  in  interaction  with  hydrogen  atoms.  A  gold  foil  within  the  poly¬ 
ethylene  is  thus  subject  to  a  flux  of  thermal  neutrons,  and  will  absorb 
them  according  to  its  cross  section  for  thermal  neutron  capture.  A 
given  flux  will  then  result  in  radioactive  gold  atoms. 
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Gold  also  has  a  resonance  cross  section  varying  from  about  1-50 
barns  for  neutrons  of  about  0.1-1  keV,  and  a  narrow  peak  of  several 
thousand  barns  at  5  eV  (Garber  and  Kinsey,  1976:400).  Since  the  flux 
integrators  thermalize  the  neutrons  to  energies  less  than  1  eV  (which  is 
below  the  resonance  interval)  and  since  the  resonance  cross  section  is 
smaller  than  the  thermal  cross  section,  the  resonance  cross  section  was 
not  considered  important  in  this  study. 

The  procedure  may  be  summarized  according  to  the  following  descrip¬ 
tion.  The  flux  integrators  are  calibrated  by  exposing  them  to  a  known 
neutron  flux,  with  gold  foils  inserted  in  the  flux  integrator.  This 
will  show  how  the  flux  integrators  respond  to  neutrons  with  a  known  flux 
and  spectrum.  Then  the  assumption  is  made  that  the  neutron  spectrum 
from  the  neutron  source  and  from  the  electron  accelerator  are  close 
enough  to  give  the  same  response  in  the  flux  integrators:  i.e.  the 
response  of  the  flux  integrators  is  not  strongly  dependent  on  the  neu¬ 
tron  spectrum.  This  is  reasonable  given  the  spectrum-flattening  charac¬ 
teristics  of  the  flux  integrators.  Their  energy  response  varies  little 
in  the  20  keV  to  20  MeV  range  (Reactor  Experiments,  Inc.,  undated:!; 
McGinley  et  al,  1976:399;  Sanders  et  al,  1984:107).  Once  calibrated, 
the  flux  integrators  are  exposed  to  the  neutron  environment  around  the 
accelerator.  From  the  foil  count  rates,  combined  with  the  counting 
efficiency  of  the  Geiger-Mueller  detector  and  the  calibration  factor  for 
the  flux  integrators,  the  neutron  flux  from  the  accelerator  can  be 
deduced. 

Huttfergg  Poiiastry 

The  important  radiation  considerations  from  neutrons  result  not 
from  primary  interaction  of  neutrons  with  tissue,  but  from  secondary 
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radiation  caused  by  the  neutrons.  Thermal  neutrons  -  that  is,  those 
below  about  0.5  eV  -  interact  chiefly  through  neutron  capture  by  atoms 
within  the  body.  These  interactions  are  capture  by  nitrogen, 
^4N(n,p)^^C,  and  capture  by  hydrogen,  ^H(n,Y)^H.  From  just  above  ther¬ 
mal  energy  to  about  10"^  MeV,  (n,p)  reactions  are  the  dominant  energy 
transfer  mechanism.  Above  this,  most  of  the  energy  and  dose  deposition 
is  by  elastic  scattering  of  hydrogen  nuclei  (Attix,  1969:468). 

HCRP  presently  recommends  a  quality  factor  (QF)  of  20  for  convert¬ 
ing  absorbed  neutron  dose  (measured  in  rads  or  centigrays)  to  equivalent 
dose  (measured  in  rems  or  centisieverts) .  Quality  factors  are  used  to 
indicate  the  expected  biological  effect  of  a  unit  of  radiation;  differ¬ 
ent  types  of  radiation  will  affect  tissue  differently  (based  on  the 
linear  energy  transfer  (LET)  of  ionizing  radiation  in  tissue). 
Alternately  if  a  factor  for  converting  neutron  fluence  to  equivalent 
dose  is  known,  the  neutron  dose  can  be  calculated  by  simply  applying 
this  conversion  factor  to  the  measured  fluence.  This  alternate  proce¬ 
dure  was  used  in  this  study. 

Meutron  Dose  Conversion 

A  number  of  studies  have  been  done  which  relate  neutron  fluence  to 
photon  dose,  neutron  dose  to  photon  dose,  and  neutron  fluence  to  neutron 
dose.  These  relations  are  usually  based  on  the  average  energy  of  the 
neutron  spectrum,  which  in  turn  depends  on  the  photon  energy  and  the 
target,  flattener,  and  collimator  assembly  materials.  For  comparison 
and  reference  a  table  of  conversion  factors  from  the  literature  is  given 
in  Appendix  A,  sorted  by  units. 
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The  Varian  Clinac  1800  runs  at  18  MeV,  Hence  it  is  reasonable  to 
expect  between  7.2x10^  n/cm^-rad^  and  1.7x10^  n/cm^-rad^  /  based  on  the 
last  two  entries  of  part  b*  of  the  table  in  Appendix  A. 

xiO^ 


£or£(MeV) 

Figure  1  Dose  Equivalent  to  Neutron  Fluence 

Figure  1  is  a  graph  of  the  relation  between  neutron  fluence  and 
neutron  dose  equivalent,  based  upon  the  energy  of  the  neutrons  (or  the 
average  energy  for  a  spectrum) •  Using  an  average  neutron  energy  of 
approximately  1  MeV  for  tungsten  photoneutrons  (Toms  and  Stephens, 
1957:77-81)  a  conversion  factor  of  2.38x10“®  remj^-n/cm2  was  used  for 
neutron  flux  to  dose  calculations  for  the  Varian  Clinac  1800  (McCall  and 
Swanson,  1979:77).  If  the  predominate  material  in  the  accelerator  head 
were  lead  (as  with  the  CGM  Saturne  1  accelerator  at  Kettering),  the 
average  photoneutron  energy  would  be  about  2  MeV  (NCRP  79,  1979:30)  and 
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a  factor  of  4* 17x10*®  remn^n/cm^  would  be  appropriate  for  the  neutron 
flux  to  neutron  dose  conversion  factor.  (NCRP  79,  1979:45;  McCall  and 
Swanson,  1979:77).  Note  these  two  conversion  factors  are  the  inverse  of 
the  factors  in  part  d.  of  the  table  in  Appendix  A.  Indeed,  part  d.  was 
derived  from  Figure  1.  These  factors  are  already  given  in  rems,  so  they 
were  used  here  as  given.  Note,  however,  that  a  set  of  quality  factors 
was  used  in  converting  D(rads)  to  H(rems).  Recent  recommendations  by 
the  ICRP  and  NCRP  are  to  double  previous  quality  factors.  If  adopted, 
all  the  doses  above  should  be  doubled. 

Once  the  neutron  fluence  is  known,  then,  the  neutron  dose  can  be 
calculated  by  the  formula 


//  =  [2. 38x10  ®rcm„Qn/cm^]<t> 


(16) 


for  the  neutrons  from  the  WPAFB  electron  accelerator  (which  uses  tung¬ 
sten).  For  measurements  done  on  the  Kettering  CGM  Saturne  1  (which  uses 
lead) ,  the  conversion  from  fluence  to  equivalent  dose  is  given  by 


//  =  [4.17  X  10‘®rem„gn/cm^)4> 


(17) 
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Ill  Method  and  Apparati 


Equipment 

The  equipment  used  for  this  study  include  the  gold  foils  and  cadmium 
covers,  the  flux  integrators,  a  Geiger-Mueller  detector  system,  radioac^ 
tive  sources,  and  the  AFIT  standard  graphite  pile.  All  this  equipment 
was  already  on  hand  at  AFIT  before  this  project  was  begun. 

Eighty-four  gold  foils  were  used.  The  supplier  assayed  them  as 
99.99%  pure.  The  diameter  and  thickness  of  the  foils  are  nominally  1 
inch  and  0.001  inch  respectively.  Their  masses  range  from  242.9  to 
300.7  mg.  Their  thicknesses  range  from  47.4  to  59.0  mg/cm^.  (Specific 
parameters  of  each  foil  are  given  in  Appendix  B« )  Tweezers  were  used  to 
handle  them  so  there  would  be  no  contamination  from  handling. 

The  cadmium  covers  are  1  inch  in  diameter,  just  large  enough  to  hold 
the  gold  foils.  They  varied  in  thickness  from  0.027  to  0.033  in.  Thit3 
thickness  is  sufficient  to  block  out  thermal  neutrons,  but  allows  most 
of  the  epithermal  (fast  or  resonant)  neutrnns  through  to  irradiate  the 
gold  foils  within  the  cadmium  co\ers*  The  effect  on  F^d  from  the  varia¬ 
tion  in  cadmium  thickness  wai  considered  to  be  less  than  other  uncer¬ 
tainties  and  therefore  neglected.  Previous  work  with  the  graphite  pile 
used  this  same  F^d  (Rossano,  1989:63;  WAOD,  1961:21),  and  these  were  the 
same  covers  used  by  Rossano  in  the  previous  study. 

The  G-M  detection  system  was  assembled  from  several  separate  pieces 
of  equipment.  At  its  core  is  a  thin  window,  pancake  style  G-M  tube. 
Model  N1002  made  by  TGM  Detectors.  It  is  filled  with  neon  and  quenched 
with  a  halogen  fill  gas.  The  window  is  mica  with  a  graphite  coating. 
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The  diameter  is  44  mm  and  the  thickness  is  about  2.0  mg/cm^*  The  mount¬ 
ing  tray  is  stainless  steel.  The  detector  and  tray  are  set  in  a  '*box" 
comprised  of  lead  bricks  to  shield  it  from  extraneous  radiation.  A 
silver  backscatterer  was  always  used  with  this  assembly. 

The  AFIT  standard  graphite  pile  consists  of  18  crossed  layers  of  4” 

X  4"  X  50"  reactor-grade  graphite.  The  pile  is  encased  in  1/32"  cadmium 
to  prevent  the  escape  or  entry  of  thermal  neutrons.  The  pile  has  9 
sliding  elements  in  which  the  gold  foils  can  be  placed  for  irradiation. 

The  neutron  source  for  the  graphite  pile,  M580,  is  a  Pu-Be  source. 
Its  output  was  1.06x10^  neutrons/second  in  September,  1989.  (The  cor¬ 
rection  to  the  output  is  shown  in  Appendix  C.  Its  output  as  specified 
by  the  manufacturer  was  8.86x10^  neutrons/ second  in  June,  1960  (Mound, 
1960:1)).  This  source  was  used  for  all  the  foil  irradiations  in  the 
graphite  pile. 

The  neurron  source  used  in  the  reactor  dome,  M1170,  is  also  a  Pu-Be 
source  with  output  of  l.lxio^  neutrons/ second  on  23  June,  1989  (calcula¬ 
tion  for  this  activity  is  in  Appendix  C) •  Its  output  as  specified  by 
the  manufacturer  was  9.00x10*  neutrons/second  on  2  March,  1962  (Mound, 
1962:1), 

The  flux  integrators  are  illustrated  in  Figure  2.  They  are  cylin¬ 
ders  with  a  diameter  of  150  mm  and  a  height  of  150  mm.  They  are  made  of 
high  density  (0.95  g/cm^)  polyethylene  with  a  0.50-mm  thick  cadmium 
cover  to  block  thermal  neutrons  from  entering  and  prevent  thermalized 
neutrons  within  from  escaping.  A  well  of  5-cro  diameter  and  7.5-cm  depth 
was  drilled  into  the  flux  integrators.  A  matching  polyethylene  plug 
with  a  cadmium  disk  at  one  end  was  made  to  fill  the  well  (Rossano, 
1989:25).  The  foils  are  inserted  into  the  bottom  of  this  well  for 
irradiations. 
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Hgtttroa  Pilft  Mgaiiurflaonti 

The  standard  graphite  pile  at  AFIT  was  used  to  calibrate  the  effi- 
cency  of  the  beta  counting  detectors.  The  graphite  pile  consists  of  18 
criss-crossed  layers  of  reactor  grade  (AGOT)  beams.  Each  beam  measures 
4**x4''x50’*  and  each  layer  is  4"  thick,  48**  wide,  and  50**  long.  The 
resulting  assembly  is  covered  on  all  sides  by  a  1/32  inch  (0.8  nun)  layer 
of  cadmium.  This  prevents  thermal  neutrons  from  escaping  the  pile  and 
also  prevents  thermal  neutrons  from  entering  the  pile.  The  cadmium 
layer  is  further  enclosed  by  an  aluminum  honey-combed  insulating  panel. 

Several  be^uns  in  the  pile  are  fashioned  to  be  easily  pulled  in  and 
out  of  the  pile.  These  beams  are  known  as  stringers.  The  neutron 
source,  M580,  fits  into  a  recess  in  the  source  stringer,  in  the  fourth 
layer  from  the  bottom  of  the  pile.  Nine  foil  stringers  ,  FS-1  through 
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FS-9,  are  located  in  every  odd  layer  beginning  from  the  bottom  layer  of 
the  pile.  Theee  stringers  have  holes  drilled  into  them  to  hold  foils; 
graphite  plugs  are  also  inserted  with  the  foils  to  space  foils  appropri¬ 
ately  and  to  surround  the  foils  with  graphite  rather  than  an  empty 
space. 

Thermal  neutron  flux  in  the  graphite  pile  is  measured  by  the  cad¬ 
mium  difference  method.  Bare  foils  are  irradiated  in  the  foil  string¬ 
ers^  and  then  cadmium-covered  foils  are  irradiated  in  the  same  location. 
The  cadmium  cover  blocks  thermal  neutrons  so  that  only  fast  neutrons 
activate  the  foils.  The  difference  in  the  activity  of  the  bare  and 
cadmium-covered  foils  is  attributed  to  thermal  neutrons,  as  was 
described  in  Chapter  II* 

Hgjpital.  MgagurcttMti 

The  flux  integrators  were  used  at  the  x-ray  cancer  therapy  machine 
at  both  the  Wright  Patterson  Medical  Center  and  at  Kettering  Memorial 
Hospital.  Measurements  were  made  both  with  and  without  patients  pres¬ 
ent*  Descriptions  are  given  in  the  following  sections* 

Measurements  Without  Patients 

A  series  of  gold  foil  irradiations  in  the  flux  integrators  was  made 
in  the  accelerator  room  with  no  patients  present.  For  the  first  set  of 
measurements,  the  top  of  the  treatment  table  was  set  at  107.5  cm  below 
the  target  of  the  machine.  This  puts  the  plane  of  the  gold  foils  at  100 
cm  from  the  target,  since  the  flux  integrators  are  15  cm  in  height  and 
the  foils  are  located  at  the  midpoint  of  the  integrator.  The  flux 
integrators  were  placed  along  the  mid-line  of  the  table,  at  distances  of 
20  cm,  30  cm,  50  cm,  and  100  cm  from  the  beam  axis.  These  were  done  two 
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at  a  time,  with  the  flux  integrators  being  set  on  opposite  sides  of  the 
beam  axis;  for  instance,  one  flux  integrator  might  be  set  with  its  cen¬ 
ter  20  cm  from  the  beam  axis  towards  the  head  of  the  table,  and  another 
set  with  its  center  100  cm  from  the  beam  axis,  but  towards  the  foot  of 
the  table.  This  prevented  one  flux  integrator  from  being  in  the  way  of 
the  other  flux  integrator  as  seen  from  the  beam  axis. 


Another  set  of  neutron  measurements  were  taken  at  four  positions 
off  of  the  table,  perpendicular  to  both  the  beam  axis  and  the  table 
length.  One  position  was  on  a  medical  stand  abutting  the  table  and  at 
the  same  height.  This  was  the  same  placement  used  for  measurements 
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during  patient  therapy.  Two  measurement  positions  were  on  chairs  2 
meters  on  either  side  of  the  table.  These  were  not  at  the  same  height 
as  the  table,  nor  were  they  at  identical  heights.  Another  position  was 
with  a  flux  integrator  set  on  a  chair  4  meters  to  one  side  of  the  table. 
Figure  3  gives  relative  locations  of  the  positions. 


Table  1:  Flux  Integrator  Locations  in  Accelerator  Room 


Position  Relative 
to 

Vertical  Mid-Plane 
of  Table 

Height  of  Plane  of 
Bottom  of  Flux 
Integrator 

Straight  Line 
Distance  from 
Target 

Next  to  table 
(on  medical  stand) 

107.5  cm  below 
plane  of  target 

106.5  cm 

2  meters  (towards 
far  wall  of  room) 

84.5  cm  above 
floor  (set  on 
table) 

232.7  cm 

2  meters  (towards 
entrance  of  room) 

63.5  cm  above 
floor  (set  on  a 
stool) 

244.13  cm 

4  meters  (towards 
entrance  of  room) 

63.5  cm  above 
floor  (set  on  a 
stool) 

423.79  cm 

The  straight  line  distances  from  the  target  head  to  these  foil 
positions  were  calculated  geometrically.  The  target  head  is  203.5  ±0.5 
cm  above  the  floor  of  the  treatment  room.  The  foil  heights  are  given  in 
Table  1.  The  difference  in  foil  height  and  target  height  defined  one 
side  of  a  triangle  for  each  position;  the  distance  from  the  center  of 
the  table  to  the  foil  positions  defined  the  base  of  such  triangles.  The 
hypotenuses  of  these  triangles  gave  the  straight-line  distances,  as 
recorded  in  Table  1.  At  the  "next  to  table"  position  the  flux  integra¬ 
tors  were  taken  as  being  centered  36.6  cm  from  the  beam  axis  (55.3 
(table  width)  +  2  -»*  7.5  (flux  integrator  radius)  +  1.5  (estimated 
distance  from  edge  of  table  to  edge  of  flux  integrator,  due  to  edge  on 
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medical  stand  which  prevented  flux  integrator  from  directly  abutting  the 
table)).  The  straight  line  distances  for  all  the  measurements  in  Table 
1  refer  to  the  distance  from  the  center  of  the  flux  integrators  (where 
the  foils  are  situated)  to  the  center  of  the  target  head,  and  may  be 
taken  as  valid  within  2  cm. 

mMiiTMienti  Witli  Patieati 

A  series  of  flux  integrator  measurements  were  made  during  cancer 
therapy  with  patients  present.  All  measurements  were  made  for  treat* 
ments  with  18  MV  x  rays.  The  accelerator  is  normally  set  at  240  monitor 
units/minute.  The  accelerator  is  calibrated  to  deliver  1  rad  per 
monitor  unit  at  depth  of  maximum  dose  in  water  (which  corresponds 
approximately  to  tissue,  since  human  tissue  is  >  80%  water).  This  depth 
is  3.2  cm  for  18  MV  x  rays. 

For  patient  treatments,  it  was  of  course  impossible  to  set  the  flux 
integrators  on  the  table  as  was  done  for  non-patient  irradiations  with 
the  accelerator.  Instead  the  flux  integrators  were  set  on  a  medical 
stand  set  at  the  same  height  as  the  table.  This  stand  was  set  as  close 
as  possible  to  the  patient  and  table,  next  to  the  area  being  irradiated. 
Here  it  was  hoped  to  capture  the  maximum  dose  of  neutrons  near  the 
patient.  (At  this  position  the  center  of  the  flux  integrator  is  36.6  ± 
2.0  cm  from  the  beam  axis,  as  calculated  in  the  previous  section.) 
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lY  ExpgriJigatflX  Romlti 


Thicknegg  Dependenciei  of  Gold  Foil  Counting 

The  thickness  of  a  foil  affects  the  activation  of  the  foil  through 
self ^shielding  and  self-absorption.  It  was  important  to  account  for 
these  processes  because  the  foils  used  in  this  procedure  were  not  infi- 
nitessimally  thin,  nor  were  they  of  exactly  uniform  thickness.  The 
majority  of  foil  activation  was  done  with  one  mil  foils,  i.e.  49  mg/cm^. 
However,  these  foils  are  only  nominally  one  mil  since  they  varied  in 
thickness  from  about  45  mg/cm^  to  about  55  mg/cm^. 

Self-shielding  occurs  during  foil  activation.  Gold  atoms  on  the 
outside  of  the  foil  block  some  neutrons  from  entering  the  foil.  This 
shields  atoms  deeper  within  the  foil  from  neutrons  that  might  activate 
them.  Self -absorption  occurs  during  counting  of  the  activated  foil. 

Beta  and  gamma  radiation  emitted  from  atoms  within  the  foil  must  travel 
through  several  layers  of  atoms  to  reach  the  surface.  Atoms  nearer  the 
foil  surface  act  to  absorb  this  outgoing  radiation;  hence  the  term  self- 
absorption.  Beta  radiation  is  absorbed  to  a  greater  extent  than  is 
gamma  radiation.  Since  the  detection  system  used  was  a  beta  counting 
system,  self-absorption  is  an  important  factor. 

It  is  rather  difficult  to  separate  self-shielding  and  self¬ 
absorption.  Experimenters  have  traditionally  combined  these  two  pro¬ 
cesses  into  a  single  correction  factor.  The  same  was  done  here. 

To  develop  this  factor,  four  half-inch  foils,  of  0.5-mil  nominal 
thickness,  were  used  in  the  standard  graphite  pile.  Four  such  foils 
will  fit  into  the  same  stringer  location  in  the  pile  without  overlapping 
each  other.  There  they  experience  the  same  neutron  flux  and  the  same 
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activation.  These  foils  were  then  counted  in  the  G-M  detector.  They 
were  counted  separately,  then  in  stacks  of  two,  three  and  then  all  four 
together.  The  saturation  count  rate  (units  of  s^^}  was  computed  by 


\C, 


(1  -e’’"*) 


(18) 


(from  Eq  (7),  i4.«C,/€).  An  expression  for  the  thickness  factor  was 

obtained  from  a  linear  regression  plot  of  the  saturation  count  rate  per 
thickness,  Cg/x  (where  x  is  thickness  in  mg/cro^)  versus  the  thickness  in 
mg/cm^  on  a  semi-log  plot.  Another  expression  for  the  thickness  factor 
was  obtained  by  plotting  the  saturation  count  rate  per  mass,  Cg/mass 
(where  mass  is  given  in  mg)  versus  the  thickness.  These  plots  are  shown 
in  Figures  4  and  S.  The  data  for  these  plots  are  in  Appendix  D.  The 
relation  for  C^/x  vs  x  was  found  to  be 

Cs 

^  -  (2.606  ±  0.036)exp[(-0.0 1  046  ±  0.00031  )a:]  (19) 


and  the  relation  for  Cg/mass  vs  x  was  found  to  be 


Cs 

- ^—-(1 .984  ±0.026)exp[(-0.0 1044  ±0.0003 1  )x]  (20) 

mass 


The  exponential  portions  of  these  terms  are  not  statistically  different 
from  the  value  found  by  a  prior  researcher.  The  value  found  previously 
is  p  «  -0.0104  1  0.0011  (Rossano,  1989:35). 
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Cs/x  [cmVnng-s] 


Figure  4  Saturation  Counts/Thickness  vs  Thickness 


Cs/Mass  [1/mg— s] 


Figure  S  Saturation  Counts/Mass  vs  Thickness 
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Datac^or  Efficiency  Dateraiinat:ion 


The  gold  foil  irradiations  in  the  graphite  pile  were  used  to  deter¬ 
mine  the  counting  efficiency  of  the  beta  counting  systems  used.  The 
detector  system  used  was  a  Geiger-Mueller  detector  with  a  silver 
backscatterer  in  the  counting  chamber.  Two  sets  of  cadmium  difference 
pairs  of  gold  foils  were  irradiated  in  stringers  1  through  6  of  the 
pile.  First  foils  were  irradiated  bare  (without  cadmium  covers)  and 
counted  in  each  detector  system  to  get  total  activity.  Then  a  similar 
set  was  irradiated  in  the  cadmium  covers  and  counted.  The  resulting 
data  gives  comparisons  between  detector  systems  and  with  previous 
results  by  Rossano. 

The  thermal  neutron  flux  at  each  stringer  location  is  well  known, 
and  thus  may  be  used  in  the  flux  equation  to  finalize  the  counting 
efficiency  term,  €,  and  any  other  constant  terms  necessary  for  the  flux 
equation  developed  earlier.  Thickness  dependent  terms  were  found  above. 
By  comparing  the  ^  ^  .erence  in  the  bare  and  cadmium  covered  foils  in 
each  stringer  ation  and  entering  all  other  known  terms  into  the  flux 
equation,  tne  only  unknovms  left  are  the  constant  terms  in  the  counting 
efficiency  term.  A  solution  for  these  constants  in  terms  of  one 
combined  constant  gives  a  gross  correction  term  to  take  care  of  any 
unforeseen  terms  that  belong  to  make  the  flux  equation  work* 

The  constant  terms  before  the  exponentials  in  Eqs  (19)  and  (20)  are 
constants  whose  value  depends  upon  what  the  foil  thickness  is  compared 
against.  The  exponential  portion  represents  the  product  of  all  the 
exponential  terms  in  the  efficiency,  c.  Thus  the  efficiency  can  be 
written  as 


e  =  /,exp(-0.01 045X) 


(21) 
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where  f  ^  is  a  constant  factor  and  the  two  exponential  values  above  have 
been  averaged. 

The  equation  for  in  the  graphite  pile,  using  cadmium  difference 
pairs,  is  given  by 


1.128/11/ c,(b)  Fc^C^jCd) 

OoNa  EffMf,  ^Cd^Cd 


(22) 


where  the  subscripts  b  and  Cd  refer  to  bare  and  cadmium  covered  foils 
respectively^  The  factor  of  1000  is  the  conversion  from  milligrams  to 
grams,  since  masses  were  measured  in  milligrams  and  AW  is  taken  as 
grams/mole.  If  a  variable  €*  -  cxp[(-0.01045)x]  is  introduced  then  from 
the  previous  paragraph  e-  .  The  equation  for  can  be  written 


I.I28/lt/1000r  C,(b)  FcdCsiCdy 

OoA^o  /x  ^Cd^  Cd 


(23) 


In  the  graphite  pile  ^ih  is  known  for  each  location.  All  other  terms  can 
be  determined  leaving  /x  as  the  only  unknown.  So  can  be  found  for 

cadmium  difference  pairs  in  each  stringer  location  in  the  pile.  Divid¬ 
ing  ^thfx  by  the  known  gives  /x  for  each  case.  Two  such  measurements 
were  done  in  each  of  stringer  locations  1  through  6  in  the  AFIT  graphite 
pile.  A  linear  least  squares  fit  was  done  of  <t>  vs  ^fx*  The  graph  of 
this  fit  is  given  in  Figure  6;  the  data  are  in  Appendix  E.  The  value  of 
f  X  from  this  fit  was  taken  as  /*  for  all  further  calculations  involving 
efficiency  c*  By  including  this  /x  it  was  found  that 

6  =  (0.428  ±  0.02 1  )exp[(-0.01 045:^:0.00031  )>:]  (24) 
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where  x  has  units  of  mg/cm^.  For  the  same  detector,  Rossano  found 
(Rossano,  1989:36) 


e  =  0.431  exp(*“0,0104x) 


(25) 


These  two  efficiency  terms  are  in  good  agreement. 

4>*fx  [neutrons/cm^ -s] 


Figure  6  Finding  fx  in  Efficiency,  e 


Increase  of  Graphite  Pile  Neutron  Cutout 

Calculation  of  the  ingrowth  of  Am-241  in  source  M580  indicates  the 
neutron  output  should  be  19.3%  higher  now  than  the  flux  recorded  in 
WADD-TR-61-174.  AF  NETF  Graphite  Standard  Pile.  The  Am-241  ingrowth  was 


0 
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calculated  according  to  a  formula  in  this  reference,  and  is  repeated  in 
Appendix  C. 

The  neutron  fluxes  for  each  stringer  location  are  given  in  Table  2, 
and  reflect  the  increased  flux  since  the  time  of  the  original  graphite 
pile  measurements.  The  original  flux  rn  1962,  as  given  in  the  WADD 
report,  is  given,  followed  by  the  expected  increase  flux  for  autumn  1989 
due  to  the  ingrowth  of  24lAm  in  source  M580.  Then  the  measured  flux  at 
each  stringer  location,  based  on  several  measurements,  is  given.  The 
nearer  the  ratio  of  the  measured  to  calculated  increase  is  to  unity,  the 
greater  our  confidence  in  the  formulas  used  for  calculating  neutron 
flux.  The  data  used  for  the  values  in  Table  2  are  given  in  Appendix  E. 


Table  2:  Neutron  Fluxes  [n/cm^-s]  in  Graphite  Pile 


Stringer 

1962 

Flux 

Increase  from 
Am-241 

Measured  Flux, 
Detector  #3 

Ratio,  1 
#3/Incr  1 

1 

3230 

3853 

3885  t  148 

1.008 

2 

6920 

8256 

7860  ±  287 

0.952 

3 

7920 

9449 

9500  ±  343 

1.005 

4 

5720 

6824 

7005  ±  251 

1.027 

5 

3180 

3794 

3905  ±  141 

1.029 

6 

1470 

1754 

1890  ±  71 

1.078 

To  ensure  the  accuracy  of  the  values  in  the  table  above,  the 
present  day  output  of  M580  was  compared  to  that  of  M1170  in  Appendix  G. 
The  original  composition  of  M1170  is  well  known  from  the  shipping  data. 
The  original  composition  of  H580  is  less  well  known  since  no  similar 
shipping  data  could  be  found  for  it  at  the  AFIT  labs.  The  likely  frac¬ 
tion  of  241pu  M580  was  derived  in  Appendix  C;  the  value  found  is 

0.0044.  If  the  correct  fractions  are  used  for  both  sources,  then  the 


31 


ratio  of  M1170  to  M580  output  as  calculated  for  a  particular  date  should 
agree  with  the  ratio  of  measured  output  for  that  date.  This  was  shown 
to  be  true  in  Appendix  G. 

rittx  Integrator  Cftlibration 

The  flux  integrators  were  calibrated  using  neutron  source  M1170. 

Its  output  during  the  time  of  experimentation  was  calculated  as  l.ixio^ 
neutrons  per  second,  as  in  Appendix  C*  The  source  and  one  or  two  flux 
integrators  at  a  time  were  hung  together  at  various  distances  (most  at 
one  meter)  from  each  other  and  at  heights  of  25,  100,  200  and  600  cm 
above  the  concrete  platform  in  the  reactor  dome  of  Building  470,  using 
the  assembly  shown  in  the  Figure  7. 

The  neutron  flux  at  the  foil  location  was  computed  by  a  l/47rr^ 
relation.  The  anisotropy  of  M1170  (1.15  at  90  degrees  orientation  (Ros- 
sano,  1989:69))  and  scatter  from  the  concrete  platform  were  also 
accounted  for.  To  account  for  scatter  from  the  platform  the  following 
equation  was  used: 


scattnrtd 
4^  dirQct 


(26) 


where  ro  is  the  straight  line  distance  from  source  to  detector  and  r^  is 
the  distance  traveled  by  a  neutron  which  undergoes  specular  reflection 
from  the  source  to  the  detector  (Jenkins,  1980:44).  Specular  reflection 
simply  means  the  angle  of  incidence  and  angle  of  reflection  are  equal. 
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E  18  the  average  neutron  energy,  in  MeV,  of  the  Bpectrum.  For  a  PuBe 
neutron  source,  the  average  energy  of  the  neutron  spectrum  is  E  «  4*2 
MeV  (Jenkins,  1980:44). 


Figure  7  Assembly  for  Source  and  Detector  in  Dome 

Since  the  neutron  flux  was  known,  and  the  saturation  activity  A,  or 

the  saturation  specific  activity  was  easily  computed  using  the  count-- 
ing  efficiency  term  found  from  the  standard  pile  data,  all  that  remained 
to  be  done  for  this  part  was  a  linear  regression  of  neutron  flux,  4>, 
versus  S.«  Such  an  equation  can  be  expressed  by 

(27) 
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where  K  represents  the  slope  and  I  is  the  intercept.  The  slope  of  this 
line,  which  is  the  calibration  constant  for  the  flux  integrators,  was 
computed  by  linear  least  squares  (Bevington,  1969:92;  Burden  and  Faires, 
1985:362-369).  The  linear  regression  and  uncertainties  were  calculated 
using  irradiation  count  numbers  302-306  and  332-335.  (The  data  for 
these  calculations  are  in  Appendix  H.)  The  resulting  equation  for  the 
flux,  given  foil  specific  saturation  activity  was  found  to  be 

({>“(2970:t  170)S.-(7.3±6.6)  (28) 

where  4> units  of  ( neutrons/ second-cro^ ) ,  S.  has  units  of  ls”^rog*^), 
and  the  constant  K  has  units  of  [neutrons-rog/cm^ ] .  (Recall  from  Eq  (14) 
that  5.  *-4. /Mass  ).  Note  the  uncertainty  for  the  intercept  is 
approximately  equal  to  the  value  for  the  intercept.  It  is  within  the 
statistical  uncertainty  of  this  value  to  allow  it  to  equal  zero.  It 
also  makes  some  intuitive  sense  here,  since  if  no  activity  is  induced 
then  one  might  expect  there  was  no  flux.  The  result  above  is  identical 
with  previous  results,  though  with  larger  uncertainties  (Rossano, 
1989:38). 
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Flux,  4>  [n/cm^-s] 


Figure  8  Least  Squares  Fit  to  Find  K(£)  for  the  Flux  Integrators 


WPAFB  Varian  Clinac  IflQQ  Measurements 

This  section  gives  results  of  a  series  of  measurements  done  at  the 
WPAFB  accelerator.  First  the  measurements  without  patients  are  pres¬ 
ented.  They  included  measurements  on  the  treatment  table  and  off  the 
treatment  table;  the  setup  for  such  measurements  was  described  in 
Chapter  III.  Following  these  are  measurements  done  during  patient  ther¬ 
apy. 

Tables  3  through  6  give  results  of  a  series  of  measurements  done 
without  patients  at  the  VfPAFB  Medical  Center.  The  Field  Size  column 
refers  to  the  size  of  the  rectangular  collimator  opening  in  the  acceler¬ 
ator  head.  Location  refers  to  distance  from  the  beam  axis.  The  dis¬ 
tance  from  the  target  to  the  table  was  107.5  cm  (so  that  the  distance 
from  the  target  to  the  foil  in  the  center  of  the  flux  integrator  was  100 
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cm).  The  flux  was  computed  by  the  formula  #  where  K  =  2970 

mg  -In/cm^  from  Eq.  (26),  and  S.  is  the  saturation  specific  activity 
[s^^mg**^]  of  the  gold  foil.  The  fluence  is  simply  the  flux  multiplied 
by  the  exposure  time*  The  dose  was  computed  using  a  conversion  factor 
of  2*38xl0~^  rem^-cm^/ neutron,  as  suggested  at  the  end  of  Chapter  II  for 
a  tungsten  target.  The  measurements  at  200  monitor  units  (M*U*)  were 
run  at  a  nominal  x-ray  machine  setting  of  240  M.U.  per  minute.  The 
remaining  measurements,  with  M.U.  readings  of  95S  to  959,  were  on  foils 
exposed  for  3  minutes  at  a  nominal  x-ray  machine  setting  of  320  monitor 
units  (M.U.)  per  minute.  The  final  readings  from  the  accelerator  con¬ 
trol  board  are  included.  Since  the  accelerator  is  calibrated  so  one 
monitor  unit  equals  one  rad^  at  depth  of  maximum  exposure  in  water,  the 
relation  between  mrem  neutron  dose  and  x-ray  absorbed  dose  was  calcu¬ 
lated  simply  by  dividing  the  dose  in  mrem^  by  the  M.U.  reading. 

Table  3  gives  measurements  done  on  the  treatment  table  with  the 
flux  integrators  centered  20  cm  from  the  beam  axis.  Field  sizes  varied 
from  10x10  cm  to  20x20  cm.  The  trend  seemed  to  show  a  larger  fluen- 
ce/radx  or  mremj^/radx  with  a  larger  field  size,  but  statistically  there 
was  no  difference  given  only  these  data  points.  The  two  measurements 
with  a  10x10  cm  field  size  averaged  (1.680±0«0S0)xl0^  n/cm^-radx  ;  the 
two  measurements  with  a  20x20  cm  field  size  averaged  (1.800±0*078)xl0^ 
n/cm^-radx  • 

Table  4  is  the  same  as  the  Table  3,  except  the  flux  integrators 
were  centered  30  cm  from  the  beam  axis,  on  the  table.  Again,  the  trend 
is  for  a  larger  fluence  or  neutron  dose  equivalent  per  x^ray  absorbed 
dose  with  larger  field  sizes,  but  statistically  there  was  no  difference 
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Table  3:  WPAFB  Accelerator  Measurements  on  Table,  20  cm 


Cnt 

# 

Field 

Size 

(cm) 

Loc’n 

(cm) 

M.U. 

Flux 

X  10-3 
n/ 6-cm2 

Fluence 

X  10-6 
n/cm2 

Dose 

(remn) 

mrem^/ 

radx 

Fluence/ 

radx 

X  10-5 

141 

10x10 

20 

200 

694 
±  55 

34.6 
±  2*7 

0*823 
±  0*064 

4.11 
±  0.32 

1*73 

±0*14 

340 

10x10 

20 

956 

865 
±  54 

155.7 
±  9*7 

3.71 
±  0.23 

3.88 

1  0.24 

1.63 

±0.10 

336 

20x16 

20 

955 

901 
±  56 

162.2 

1  10.8 

3.86 
±  0.26 

4.04 
±  0.27 

1.70 
±  0.11 

342 

20x20 

20 

959 

959 
±  60 

172.6 
±  10.8 

4*11 
±  0*26 

4*29 
±  0.27 

1*80 
±  0.11 

346 

20x20 

20 

956 

956 

1  60 

172.1 
±  10*8 

4.10 
±  0.26 

4*29 
±  0*27 

1.80 
±  0*11 

with  different  field  sizes  from  this  data*  The  average  with  a  10x10  cm 
field  size  was  (1.472±0*025)xl0^  n/cm^-radx  ;  with  a  20x20  cm  field  size 
the  average  was  (1* 567±0.025)xl0^  n/cm^*radx  • 


Table  4:  WPAFB  Accelerator  Measurements  on  Table,  30  cm 


Cnt 

# 

Field 

Size 

[cm] 

Loc’ 

(cm) 

M.U* 

Flux 

X  10-3 
n/s-cro^ 

Fluence 

X  10-6 
n/cttfi 

Dose 

(remn) 

mremf)/ 

radx 

Fluence/ 

radx 

X  10“5 

145 

10x10 

30 

200 

611 
±  41 

30*4 
±  2*0 

0*724 
±  0*048 

3.62 
±  0.24 

1.52 
±  0.10 

146 

10x10 

30 

200 

576 
±  38 

28*7 
±  1.9 

0.683 
±  0.045 

3.42 
±  0.22 

1.435 
±  0.095 

338 

10x10 

30 

958 

778 
±  48 

140.0 
±  8.6 

3*33 
±  0*20 

3.48 
±  0*21 

1.461 
±  0.090 

344 

20x20 

30 

958 

847 
±  53 

152*5 
±  9.5 

3*63 
±  0*22 

3*79 
±  0.23 

1*592 
±  0*099 

348 

20x20 

30 

957 

820 
±  51 

147*6 
±  9.2 

3*51 
±  0*26 

3*67 
±  0*27 

1.542 
±  0.096 
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Table  5  gives  measurements  on  the  table  when  the  flux  integrators 


were  centered  50  cm  from  the  beam  axis  on  the  table*  The  trend  of 
larger  fluence/radx  or  mrem^/radx  with  larger  field  sizes  reverses  at 
this  distance;  the  neutron  fluence  and  neutron  eG[uivalent  dose  per  x-ray 
absorbed  dose  decreases  with  larger  field  sizes*  However,  again  there 
is  no  real  difference  statistically*  At  10x10  cm  field  size  the  average 
is  (1* 304310 *0083)xl0^  n/cm^-radx  $  and  at  20x20  cm  field  size  the  aver¬ 
age  is  (1*2765±0*0075)x105  n/cm^-radx  • 


Table  5:  WPAFB  Accelerator  -  Measurements  on  Table,  50  cm 


Cnt 

# 

Field 

Size 

(cm) 

LoC 

(cm] 

M*U* 

Flux 

X  10-3 

n/B-cm2 

Fluence 

X  10“6 
n/cm2 

Dose 

(remnl 

mrem^/ 

radx 

Fluence/ 

radx 

X  10“5 

142 

10x10 

50 

200 

528 
±  35 

26.3 
i  1.8 

0*626 
±  0*043 

3*13 
±  0*22 

1.315 
±  0.090 

144 

10x10 

50 

200 

511 
±  35 

26*2 
±  1*8 

0*624 
±  0*043 

3*12 
±  0*22 

1*310 
±  0.090 

341 

10x10 

50 

956 

684 
±  43 

123*1 
±  7*7 

2*93 
±  0*18 

3*06 
±  0*19 

1*288 
±  0.081 

337 

20x16 

50 

955 

672 
±  42 

HPRWQIII 

110019 

2*88 
±  0*18 

3.02 
t  0.19 

1*267 
±  0*080 

343 

20x20 

50 

959 

684 
±  43 

123*1 
±  7*7 

2.93 
±  0.18 

3*06 
±  0.19 

1*284 
±  0*080 

347 

20x20 

50 

959 

674 
±  42 

2.88 
±  0.18 

3*01 
±  0*19 

1*269 
±  0*079 

The  final  table  in  this  series.  Table  6,  gives  results  for  on-table 
measurements  when  the  flux  integrators  were  centered  100  cm  from  the 
beam  axis*  The  fluence  per  radx  and  mremi^  per  radx  again  show  decreas¬ 
ing  trends  with  larger  field  sizes,  but  statistically  there  is  not  a 
difference*  At  10x10  cm  field  size  the  average  was  (0*9315±0*0085)xl0^ 
n/cm^-radx  ;  at  20x20  cm  field  size  the  average  was  (0*902510*0065)xl0^ 
n/cm^-radx  • 
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Table  6:  WPAFB  Accelerator  *  Measurements  on  Table,  100  cm 


Cnt 

# 

Field 

Size 

(cm) 

LoC 

(cm) 

M.U. 

Flux 

X  10-3 
n/s-cm^ 

Fluence 

X  10-6 
n/cm2 

Dose 

(reran) 

mremn/ 

radx 

Fluence/ 

radx 

X  10-5 

147 

10x10 

100 

200 

381.3 
±  26.6 

19.0 
±  1.3 

0.452 
±  0.031 

2.26 
±  0.16 

0.940 
±  0.065 

339 

10x10 

100 

958 

491 
±  31 

88.4 
±  5.6 

2.10 
±  0.13 

2.19 
±  0.14 

0.923 
±  0.058 

345 

20x20 

100 

958 

484 
±  30 

2.07 
±  0.13 

2.16 
±  0.14 

0.909 

1  0.056 

349 

20x20 

100 

957 

476 
±  30 

2.04 
±  0.13 

2.13 
±  0.14 

0.896 
±  0.056 

These  four  tables  (3*6)  show  the  neutron  fluence  per  photon  rad  to 
be  within  the  range  established  by  previous  literature.  Specifically, 
these  measurements  on  the  treatment  table  range  from  a  high  of  1.80x10^ 
n/cm2-radx  at  20  cm  axial  distance  from  the  beam  to  a  low  of  8.96x10^ 
n/cro2-rad^  at  100  cm  axial  distance  from  the  beam,  both  at  20x20  cm 
field  size.  A  range  of  7.4x10^  n/cm^-rad^  to  1.7x10^  n/cm^-rad^  was 
expected  from  previous  experiment,  as  from  Table  16.  There  may  be  a 
trend  between  neutron  measurements  and  field  size.  At  closer  distances 
there  seem  to  be  more  neutrons  released  when  the  field  size  is 
increased,  and  at  longer  distances  there  seem  to  be  fewer  neutrons  with 
a  larger  field  size. 

The  next  four  tables  (7-10)  include  more  irradiations  at  the  Varian 
Clinac  1800,  again  without  patients.  In  these  the  flux  integrators  were 
not  set  on  the  treatment  table.  In  the  first  two  tables  they  were  set 
on  the  medical  stand  next  to  the  table,  to  simulate  their  placement 
during  measurements  with  patients.  Recall  that  this  location  is  cen¬ 
tered  at  36.6  cm  from  the  beam  axis  (Chapter  III). 
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From  Table  7  the  neutron  fluence  averaged  ( 1. 149±0.018)xl0^ 
n/cm^-rad^  at  the  aide-of-table  position  for  a  10x10  cm  field  size.  For 
a  20x20  cm  field  size  the  fluence  averaged  (1.235±0.010)xloS  n/cm^-radx* 
In  this  case  there  is  a  statistical  difference  between  the  neutron  flu¬ 
ence  per  radx  for  the  different  field  sizes — the  neutron  fluence  is 
higher  for  a  larger  field  size. 


Table  7:  WPAFB  Accelerator  -  Measurements  Beside  Table 


Cnt 

# 

Field 

Size 

(cmj 

M.U. 

Flux 

X  10*3 

n/e-cjn2 

Fluence 

X  10-6 
n/cm2 

Dose 

(mremn) 

mremjj 

radx 

Fluence/ 

radx 

X  10“5 

310 

10x10 

317 

601 
±  39 

36.1 
±  2.3 

859 

1  55 

2.71 
±  0.17 

1.139 
±  0.073 

311 

10x10 

316 

624 
±  39 

37.4 
±  2.3 

890 
±  55 

2.82 
±  0.17 

1.184 
±  0.073 

312 

10x10 

316 

653 
±  41 

39.2 

1  2.5 

933 
t  56 

2.95 
t  0.18 

1.124 
±  0.079 

313 

20x20 

318 

667 

1  44 

952 

1  62 

2.99 
±  0.19 

1.258 
t  0.082 

314 

20x20 

313 

656 
±  41 

39.4 
±  2.5 

938 
±  56 

3.00 
±  0.18 

1.259 
±  0.080 

315 

20x20 

318 

647 
i  41 

38.8 
t  2.5 

923 
±  56 

2.90 
i  0.18 

1.220 
±  0.079 

320 

20x20 

1914 

652 
±  40 

235 
t  14 

5590 
±  330 

2.92 
±  0.17 

1.228 
±  0.073 

324 

20x20 

1917 

645 
±  40 

232 
t  14 

5520 
±  330 

2.88 
t  0.17 

1.210 
±  0.073 

Table  8  is  the  same  as  Table  7,  except  a  steel  wedge  and  a  lead 
block  were  used  to  shape  the  beam,  as  is  often  done  in  treatment*  Count 
numbers  316  and  317  were  done  with  a  45^  steel  wedge  in  the  gantry  head. 
Count  numbers  318,  319,  and  328  were  done  with  a  half-beam  lead  block 
set  at  0^  orientation  in  the  gantry  head.  These  data  average 
(1. 16410. 025 )xl0^  n/cm^-radx  for  the  steel  wedge,  and  ( 1.26510.018 )xl0^ 


40 


n/cm^-radx  for  the  lead  block.  There  la  no  etatiatical  difference 
between  each  of  these  averages  and  the  average  for  the  same  condition 
without  a  wedge  or  block  (Table  7),  though  the  trend  seems  to  be  fewer 
neutrons  with  the  steel  wedge  and  more  neutrons  with  the  lead  block. 
However  there  is  a  statistical  difference  between  using  the  block  and 
using  the  wedge;  the  difference  in  the  averages  is  (0. 08610. 021)xl0^ 
n/cm2-radx  . 


TaUble  8:  WPAFB  Accelerator  -  Measurements  Beside  Table  with  Steel  or 
Lead  Beam  Shapers;  20x20  cm  Field  Size 


Cnt 

# 

Beam 

Shaper 

M.U. 

Flux 
x  10*-3 
n/s-cm2 

Fluence 

X  10-6 
n/cm2 

Dose 

tnireron] 

mrenin 

/ 

radx 

Fluence/ 

radx_ 

X  10-5 

316 

45”  Btsel 
wedge 

317 

602 
±  40 

859 
±  57 

2.71 
±  0*18 

1*139 
±  0*076 

317 

45^  steel 
wedge 

318 

630 
±  40 

37*8 
±  2*4 

900 
±  57 

2.83 
±  0.18 

1.189 
±  0.075 

318 

Pb  block 

316 

677 
±  44 

40*6 
±  2*6 

966 
±  62 

3*06 
±  0.20 

1.285 
±  0*082 

319 

Pb  block 

318 

679 
±  43 

40.7 
±  2.6 

969 
±  62 

3.07 
±  0.19 

1*280 
±  0*082 

328 

Pb  block 

1913 

654 
±  41 

235 
t  15 

5590 
±  360 

2.92 
±  0*19 

1.229 
t  0.078 

In  TaUDle  9  the  flux  integrators  were  set  along  a  line  perpendicular 
to  both  the  axis  of  the  table  and  the  axis  of  the  x-ray  beam*  These 
were  at  2  and  4  meters  from  the  center  of  the  table.  The  measurements 
located  at  **2,  wall"  were  set  at  2  meters  towards  the  far  wall  of  the 
treatment  room,  as  described  in  Table  2  in  chapter  3*  Those  located  at 
"2,  door"  and  "4,  door"  were  set  2  meters  and  4  meters,  respectively, 
towards  the  door  of  the  treatment  room*  Additionally  count  numbers 
329-331  were  done  with  the  lead  block  in  the  gantry  head,  as  mentioned 
in  the  previous  table.  No  statistical  difference  in  the  measurements 
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could  be  seen  from  these  data  with  the  inclusion  of  the  lead  block. 

There  were  more  neutrons  measured  towards  the  wall  of  the  treatment  room 
than  towards  the  door;  one  would  expect  more  neutron  reflection  on  that 
side  since  the  corner  of  the  room  and  the  walls  are  closer  on  that  side. 


Table  9:  WPAFB  Accelerator  -  Measurements  Off  Table 


Cnt 

# 

Field 

Size 

[cm] 

Loc'n 

(m] 

M.U. 

Flux 

X  10“3 
n/B-cm2 

Fluence 

X  10-® 
n/cin^ 

Dose 

(mremn) 

mrem^ 

/ 

radx 

Fluence/ 

X  10-4 

321 

20x20 

2, 

wall 

1914 

270 
±  17 

97.2 
±  6.1 

2310 
±  150 

1.207 
±  0.078 

5.08 
±  0.32 

325 

20x20 

2, 

wall 

1917 

278 
±  18 

2380 
±  150 

1.242 
±  0.078 

5.22 
±  0.34 

329 

20x20 

2, 

wall 

1913 

290 
±  18 

2480 
±  150 

1.296 
±  0.078 

5.46 
±  0.34 

322 

20x20 

2, 

door 

1914 

220 
±  14 

79.2 
±  5.0 

1880 
±  120 

0.982 
i  0.063 

4.14 
±  0.26 

326 

20x20 

2, 

door 

1917 

227 
±  14 

1940 
±  120 

1.012 
±  0.063 

4.26 
±  0.26 

330 

20x20 

2, 

door 

1913 

228 
±  14 

1950 
±  120 

1.019 
±  0.063 

4.29 
±  0.26 

323 

20x20 

4, 

door 

1914 

109.0 
±  7.2 

39.2 
±  2.5 

933 
±  60 

0.487 
±  0.031 

2.05 
±  0.13 

327 

20x20 

4, 

door 

1917 

95.3 
±  6.3 

34.3 
±  2.3 

816 
±  55 

0.426 

1  0.029 

1.79 
±  0.12 

331 

20x20 

4, 

door 

1913 

99.6 
±  6.5 

35.6 
±  2.3 

847 
±  55 

0.443 
±  0.029 

1.86 

1  0.12 

Table  10,  the  final  table  in  this  section,  presents  results  of  flux 
integrator  measurements  with  patients  during  their  cancer  therapy.  The 
conditions  were  necessarily  less  well  controlled  for  these  measurements 
than  for  those  without  patients.  The  field  sizes  were  not  constant;  a 
variety  of  lead  blocks  were  used  to  shape  the  x^-ray  beam  for  individual 
treatments;  the  exposure  times  were  short  in  comparison  to  the  measure¬ 
ments  without  patients;  the  position  of  the  flux  integrators  was  less 


42 


well  controlled.  For  instance,  count  #298  was  placed  by  the  patient's 
calf  rather  than  beside  the  treated  area  because  of  treatment  conditions 
(the  gantry  had  to  rotate  around  the  table  and  the  flux  integrator  had 
to  be  kept  out  of  the  way);  this  gave  a  lower  neutron  dose  rate  for  this 
measurement,  since  the  flux  integrator  was  further  away  than  in  other 
measurements.  For  the  other  measurements  the  neutron  equivalent  dose 
per  x-ray  absorbed  dose  (mremn/rad^)  results  vary  somewhat,  but  is  of 
the  seune  order  of  magnitude  as  the  without -patients  measurements.  This 
demonstrates  the  validity  of  this  technique  for  measuring  neutron  dose 
for  x-ray  treatments  at  the  WPAFB  facility. 


Table  10:  WPAFB  Accelerator  -  Measurements  with  Patients 


Cnt 

# 

Field 

Size 

I  cm) 

Monitor 

Units 

Flux 

X  10'3 

n/s-cm^ 

Fluence 

X  10-6 

n/cin2 

Dose 

(mremn) 

mremn 

radx 

149 

294 

852  ±  55 

47.0  ±  3.0 

1119  ±  71 

3.81  ±  0.24 

150 

200 

470  ±  33 

18.0  ±  1.3 

428  ±  31 

2.14  ±  0.16 

151 

238 

679  t  45 

30.1  1  2.0 

716  ±  48 

3.01  ±  0.20 

152 

318 

924  ±  60 

44.4  ±  2.9 

1057  ±  69 

3.32  t  0.22 

298 

8x8 

243 

26.1  ± 
3,9 

2.79  i  0,42 

66.4  t  1.0 

0.2733  ± 
0.0041 

299 

14x25, 

15x25 

238 

641  ±  40 

28.8  ±  1.8 

685  ±  43 

2.88  ±  0.18 

307 

10x17, 

12x17 

199 

682  ±  43 

25.8  ±  1.6 

614  1  38 

3.09  ±  0.19 

300 

17x17, 

14x17 

258 

358  ±  23 

17.6  ±  1.1 

419  ±  26 

1.62  ±  0.10 

308 

17x17, 

14x17 

258 

380  ±  24 

18.9  ±  1.2 

450  ±  29 

1.74  ±  0.11 

301 

17x20, 

16x20 

220 

607  ±  39 

24.8  1  1.6 

590  ±  38 

2.68  ±  0.17 

309 

17x20, 

16x20 

220 

463  ±  29 

19.4  ±  1.2 

462  ±  29 

2.10  ±  0.13 
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From  Table  10  it  is  seen  the  neutron  equivalent  dose  per  x-ray 
absorbed  dose  varied  from  a  low  of  1.62  mrem^/radx  to  a  high  of  3.81 
mremn/radx  #  or  dose  to  patient  of  0.42  to  1.12  remn  • 

Kgttoriag-CflM-Sfttttrag  1  Mcaguroagnti 

A  series  of  measurements  were  also  made  at  Kettering  Memorial  Hos¬ 
pital.  The  Kettering  x-ray  machine  is  a  French-made  CGM  Saturne  1, 
which  operates  at  10  MV*  It  is  magnetron  driven,  with  a  beeun  energy  of 
14-15  MeV*  It  has  a  lead  collimator,  so  a  flux-to-dose  conversion  fac¬ 
tor  of  4*17x10^^  remn-cm^ /neutron  was  used.  The  setup  was  essentially 
identical  to  that  done  at  the  WPAFB  Medical  Center. 

Tables  11-14,  following,  give  the  results  of  a  series  of  measure¬ 
ments  at  Kettering  Memorial  Hospital.  The  tables  give  measurements  at 
20,  30,  50,  and  100  cm  on  the  treatment  table,  along  its  axis.  The  CMC 
Saturne  1  was  run  at  1500  monitor  units,  so  the  doses  were  all  divided 
by  1500  to  obtain  the  relation  between  mremn  and  radx* 

Table  11  is  for  the  measurements  at  20  cm.  Two  field  sizes,  10x10 
cm  and  20x20  cm,  were  used.  The  measurements  of  fluence/radx  and 
mremn/radx  were  statistically  larger  for  the  larger  field  size  at  this 
distance:  averages  were  (5. Il±0.21)xl0^  n/cm^-radx  for  10x10  cm  field 
size  vs.  (5.98S±0.03S)xl0^  n/cm^-radx  ^or  20x20  cm  field  size.  It  seems 
neutrons  may  be  blocked  at  closer  distances  by  a  smaller  field  size. 

Table  12  gives  measurements  at  30  cm  axial  distance  on  the  treat¬ 
ment  table.  Again,  measurements  were  made  with  field  sizes  at  both 
10x10  cm  and  20x20  cm.  The  larger  field  size  again  gives  larger  values 
for  the  neutron  measurements  per  x  ray*  Averages  were  (4.72±0.06)xl04 
n/cm2-radx  st  10x10  cm  and  (S.18±0.16)xl0^  n/cm^-radx  at  20x20  cm. 
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Table  11:  Kettering  CGM  Saturne  1  -  Measurements  on  Table,  20  cm 


Cnt 

# 

Field 

Size 

(cm] 

Loc '  n 
(cm] 

Flux 

X  10-3 
n/B-cm2 

Fluence 

X  10-6 
n/cm2 

Dose 

(reran) 

rareran/ 

radx 

Fluence/ 

radx 

X 

350 

10x10 

20 

198.4 
t  3.6 

73.6 
±  1.3 

3.069 
±  0.054 

2.046 
±  0.036 

4.907 
±  0.087 

355 

10x10 

20 

mgm 

79.8 
±  1.4 

3.328 
±  0.058 

2.219 
±  0.039 

5.320 
±  0.093 

362 

20x20 

20 

243.5 
±  4.8 

90.3 
±  1.8 

3.766 
±  0.075 

2.511 
±  0.050 

6.02 
±  0.12 

370 

20x20 

20 

89.3 
±  1.7 

3.724 
±  0.071 

2.483 
±  0.047 

5.95 
t  0.11 

Table  12:  Kettering  CGM  Saturne  1  -  Measurements  on  Table,  30  cm 


Cnt 

# 

Field 

Size 

(era) 

Loc’n 

(era) 

Flux 

X  10-3 

n/B-cra2 

Fluence 

X  10-6 
n/cra2 

Dose 

(reran) 

rareran/ 

radx 

Fluence/ 

radx 

X  10-4 

353 

10x10 

30 

■m 

2.990 
±  0.054 

1.993 
±  0.036 

4.780 
±  0.087 

359 

10x10 

30 

188.6 
±  3.6 

69.9 
±  1.3 

2.915 
±  0.054 

1.943 
±  0.037 

4.660 
±  0.087 

366 

20x20 

30 

mssBrn 

MBM 

3.344 
±  0.067 

2.229 
±  0.045 

5.347 
±  0.011 

374 

20x20 

30 

75.3 
±  1.6 

3.140 
±  0.067 

2.093 
±  0.045 

5.020 
±  0.011 

Table  13  shows  measurements  when  the  flux  integrators  were  centered 
50  cm  from  the  beam  axis  on  the  table.  Both  10x10  cm  and  20x20  cm  field 
sizes  were  used.  At  this  distance,  however,  there  was  no  statistical 
difference  in  the  neutron  measurements  per  rad  x-ray  with  varying  field 
size.  Perhaps  neutron  scatter  and  reflection  has  begun  to  become  more 
uniform  with  larger  distances  from  the  target. 


45 


Table  13:  Kettering  CGM  Saturne  1  -  Measurements  on  Table,  50  cm 


o 

ft 

Field 

Size 

(cmj 

Loc  'n 
(cm] 

Flux 

X  10-3 
n/s-cm^ 

Fluence 

X  10-6 
n/cm2 

Dose 

(remnJ 

mreroji/ 

radx 

Fluence/ 

radjt 

X  10-4 

351 

10x10 

50 

59.9 
±  1.2 

2.498 
t  0.050 

1.665 
±  0.033 

3.993 
±  0.080 

358 

10x10 

50 

MBHI 

59.9 
±  1.2 

2.498 
±  0.050 

1.665 
±  0.033 

3.993 

1  0.080 

363 

20x20 

50 

^BBDI 

2.577 
±  0.054 

1.718 
±  0.036 

4.120 
t  0.087 

371 

20x20 

50 

156.2 
±  42 

2.414 
t  0.050 

1.609 
±  0.033 

3.860 
±  0.080 

Table  14  gives  the  on-table  measurements  when  the  flux  integrator 
was  centered  100  cm  from  the  beam  axis.  Here  the  trend  seems  to  agree 
with  the  WPAFB  measurements,  in  that  the  neutron  fluence  and  equivalent 
dose  per  x-ray  absorbed  dose  decreased  with  a  larger  field  size  at  the 
further  distances.  The  statistical  difference  is  not  strong,  however, 
based  on  this  data.  The  averages  were  (2.99±0.14)xl0^  n/cra^-rad^  and 
(2 . 720±0.048)xl0^  n/cm^-rad^  at  10x10  cm  and  20x20  cm  respectively. 

For  a  10  MV  x-ray  machine  a  neutron  fluence  to  photon  dose  measure¬ 
ment  on  the  order  of  1x10^  n/cm^-rad^  would  be  expected  from  Table  16. 
Table  11  shows  averages  of  (5. Il±0.21)xl0^  n/cm^-rad^  for  a  10x10  cm 
field  and  (5.985±0.035)xl04  n/cm^-radx  for  a  20x20  cm  field,  both  at  20 
cm  distance  from  the  beaun.  These  are  perhaps  3^  times  larger  than 
expected,  but  only  a  limited  number  of  references  for  10  MV  accelerators 
were  included.  Without  more  knowledge  of  typical  10  MV  accelerators  it 
cannot  be  said  here  whether  these  are  truly  larger  measurements  than  the 
norm.  This  accelerator  does  produce  a  lower  fluence  per  x-ray  dose  than 
the  WPAFB  accelerator. 
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Taible  14:  Kettering  CGK  Saturne  1  -  Measuremente  on  Table,  100  cm 


Cnt 

# 

Field 

Size 

[cm] 

Loc'n 

[cm] 

Flux 

X  10-3 

n/B-cm2 

Fluence 

X  10-6 
n/cm2 

Dose 

[remn] 

mremji/ 

radx 

Fluence/ 

radx 

X  10-4 

352 

10x10 

100 

115.3 
t  2.6 

42.75 
±  0.96 

1.728 
±  0.040 

1.152 
±  0.027 

2.850 
±  0.064 

354 

10x10 

100 

126.3 
±  2.7 

46.8 
±  1.0 

1.952 
±  0.042 

1.301 
±  0.028 

3.120 
t  0.067 

367 

20x20 

100 

111.9 
±  2.8 

41.5 
±  1.0 

1.731 
±  0.042 

1.154 
±  0.028 

2.767 
±  0.067 

375 

20x20 

100 

108.1 

40.08 
±  0.93 

1.671 
±  0.039 

1.114 
±  0.026 

2.672 
t  0.062 

Table  15  includes  more  measurements  done  at  the  Kettering  accelera¬ 
tor.  Measurements  at  200  and  400  cm  distance  from  the  beam  axis  were 
done  perpendicularly  to  the  table  axis,  with  the  flux  integrators  placed 
on  chairs.  The  measurements  at  2  meters  were  done  with  the  bottom  of 
the  flux  integrator  at  62.5  cm  above  the  floor.  At  4  meters,  count 
number  356  was  done  with  the  bottom  of  the  flux  integrator  at  61.5  cm 
from  the  floor;  count  number  357  was  done  with  the  bottom  of  the  flux 
integrator  at  60.5  cm  above  the  floor.  All  other  measurements  at  4 
meters  were  done  with  the  bottom  of  the  flux  integrators  at  62.5  cm 
above  the  floor.  Unfortunately  the  height  of  the  accelerator  head  tar¬ 
get  is  unknown  to  this  researcher,  so  the  straight  line  distance  from 
the  target  to  the  foils  within  the  flux  integrators  could  not  be 
calculated.  These  measurements  are  included  for  comparison  only.  There 
seems  to  be  a  trend  of  higher  neutron  measurements  at  smaller  field 
sizes,  consistent  with  the  trends  seen  at  50  and  100  cm  on  the  table. 
However  no  statistical  difference  between  the  field  sizes  could  be 
shown.  Averages  at  200  cm  were  (1.344±0.049)xl0^  n/cm^-rad^  at  10x10  cm 
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Table  15:  Kettering  COM  Saturne  1  -  Measurements  Off  Table 


ft 

Field 

Size 

[cm] 

Loc'n 

tcm) 

Flux 

X  10“3 
n/e-cm^ 

Fluence 

X  10-6 
n/cm2 

Dose 
[ reran] 

mremn/ 

radx 

Fluence/ 

radx 

X  10-4 

356 

10x10 

200 

56.3 

1  1.6 

20.90 
±  0.S9 

0.872 
±  0.025 

0.581 
±  0.017 

1.393 
±  0.039 

360 

10x10 

200 

52.4 

1  1.5 

19.43 

1  0.56 

0.810 
±  0.023 

0.540 
±  0.015 

1.295 
±  0.037 

364 

20x20 

200 

49.9 

1  1.5 

18.50 
±  0.56 

0.771 
±  0.023 

0.514 
±  0.015 

1.233 
±  0.037 

368 

20x20 

200 

18.87 
±  0.56 

0.787 
±  0.023 

0.525 
±  0.015 

1.259 
±  0.037 

372 

20x20 

200 

18.69 
±  0.56 

0.779 
±  0.023 

0.519 
±  0.015 

1.246 
±  0.037 

376 

20x20 

200 

48.7 

1  1.8 

18.06 
±  0.67 

0.753 
±  0.028 

0.502 
±  0.019 

1.204 
±  0.045 

357 

10x10 

400 

23.6 

1  1.1 

58.75 
±  0.41 

0.365 
±  0.017 

0.243 
±  0.011 

5.833 
±  0.027 

361 

10x10 

400 

25.6 

1  1.2 

9.49 
±  0.44 

0.396 
±  0.018 

0.264 
±  0.012 

6.327 
±  0.029 

365 

20x20 

400 

BBQI 

8.01 
±  0.41 

0.334 
±  0.017 

0.223 
±  0.011 

5.340 
±  0.027 

369 

20x20 

400 

23.2 

1  1.1 

8.60 
±  0.41 

0.359 
±  0.017 

0.239 
±  0.011 

5.733 
±  0.027 

373 

20x20 

400 

mm 

9.08 
±  0.44 

0.379 
±  0.018 

0.253 
±  0.012 

6.053 
±  0.029 

377 

20x20 

400 

22.8 
±  1.2 

8.45 
±  0.44 

0.352 
±  0.018 

0.235 
±  0.012 

5.633 

1  0.029 

and  (1*23510. 012 )xl04  n/cm^-rad^  at  20x20  cm.  At  400  cm  the  averages 
were  (6.08±0.25)xl04  n/cm^-radx  at  10x10  cm  and  (5.6910. 15)xl04  n/cm^- 
radx  at  20x20  cm. 
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V  ConcliiBlQii 


The  first  step  in  this  project  was  using  the  standard  graphite  pile 
to  ascertain  the  counting  efficiency  of  the  Geiger-Mueller  detector  used 
to  count  activated  gold  foils.  Foils  with  a  nominal  thickness  of  0.5 
mil  were  irradiated  and  counted  individually  and  in  stacks  of  two  and 
three  (to  make  nonimal  thicknesses  of  1  mil  and  1.5  mil)  to  give  a  graph 
of  counts  vs.  thickness.  The  resulting  exponential  relation  was 
combined  with  a  series  of  cadmium-*dif ference  measurements  to  obtain  the 
relation  for  efficiency#  e,  to  use  when  counting  foil  activity  with  the 
detector.  This  relation  was  found  to  be 


€  =  (0.428*0.021  )exp[(- 0.01 045  ±0.0003  l)x]  (24) 

where  x  has  units  of  mg/cm^  .  The  efficiency  term  c  is  unitless. 

The  next  major  step  was  calibration  of  the  flux  integrators.  A 
series  of  measurements  at  various  heights  above  the  concrete  platform  in 
the  AFIT  reactor  dome  gave  the  relation 

(1>*(2970*  180)S^  (28a) 


which  was  used  for  all  subsequent  irradiations  made  with  gold  foils  in 
the  flux  integrators.  5.  has  units  of  It  is  derived  from 

A^/Mass  where 
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A 


(10) 


\c 


nat 


Xr 

e 


€(1 


and  te#  tc#  and  ty,  are  the  exposure  time,  count  time,  and  wait  time  of 
the  gold  foil,  respectively «  To  complete  the  summary  of  important  equa¬ 
tions  the  relation  is  repeated,  where 


(1 


(18) 


Both  i4.  and  Cg  have  units  of 

The  final  part  of  this  thesis  study  involved  measurements  done  with 
the  flux  integrators  without  and  with  patients  present  at  the  WPAFB 
Varian  Clinac  1800  x-ray  therapy  machine*  Measurements  were  also  done 
at  the  Kettering  CGM  Saturne  1  10  MV  cancer  therapy  machine  for  compari¬ 
son. 

Results  and  Further  Suggestions 

In  the  saturation  activity  equation  and  saturation-to-f lux  equation 
developed  here,  the  constants  for  detector  efficiency  e  and  flux  inte¬ 
grator  calibration  constant  R  were  not  significantly  different  frcxn 
those  determined  in  the  previous  work  by  Rossano.  This  confirms  the 
stability  of  the  measuring  system  developed  by  him  and  gives  greater 
confidence  of  its  validity  and  usefulness. 

Application  of  the  procedure  at  the  Varian  Clinac  1800  at  the  WPAFB 
Medical  Center  gave  fluences  within  the  range  of  previously  published 
results,  as  in  Table  16.  Application  at  the  CGM  Saturne  1  at  Kettering 
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Memorial  Hoapltal  further  proved  the  usability  of  this  scheme*  The 
procedure  outlined  here  should  be  quite  feasible  for  periodic  monitoring 
of  the  neutron  flux  of  the  accelerator  at  the  VTPAFB  Medical  Center* 

Flux  integrator  irradiations  at  the  WPAFB  Medical  Center  without 
patients  gave  (1.80±0*ll)xl0^  n/cm^^-radx  for  measurements  on  the  treat¬ 
ment  table,  20  cm  from  the  beam  axis,  with  a  field  size  of  20x20  cm* 

The  corresponding  average  neutron  equivalent  dose  per  x-ray  absorbed 
dose  was  4 *29010. 010  mremn/radx  •  Measurements  at  a  position  immedi¬ 
ately  beside  the  treatment  table  averaged  (1.235±0.035)xl0^  n/cm^-radx 
with  a  field  size  of  20x20  cm;  the  corresponding  average  neutron 
equivalent  dose  per  x-ray  absorbed  dose  was  2*93610*080  mremn/radx  • 
Measurements  were  also  done  at  30,  SO,  and  100  cm  from  the  beam  axis  on 
the  treatment  table,  as  well  as  2  meters  from  the  beam  axis  on  either 
side  of  the  table  and  at  4  meters  from  the  beam  axis  towards  the  door  of 
the  treatment  room. 

Based  on  the  measurements  at  different  locations  with  different 
field  sizes,  there  seems  to  be  an  interesting  trend*  For  measurements 
done  closer  than  about  40  cm  from  the  beam  axis,  there  is  a  larger 
neutron  fluence  with  larger  field  size*  For  measurements  at  50  cm  and 
further,  there  is  a  larger  neutron  fluence  with  smaller  field  size.  A 
smaller  field  size  means  there  is  more  material  in  the  path  of  the  beam; 
we  might  expect  more  neutrons  to  be  released*  However,  this  same  mate¬ 
rial  seems  to  deflect  neutrons  from  the  area  close  to  it*  This  could 
explain  why  there  are  fewer  neutrons  measured  close  to  the  beam  and  more 
further  away  when  the  field  size  is  reduced* 

The  difference  in  fluence/radx  between  the  steel  wedge  and  the  lead 
block  (Table  8)  suggests  it  might  be  safer  to  the  patient  to  use  steel 
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for  beam  shaping.  However,  a  variety  of  factors  which  surely  contrib¬ 
uted  to  the  difference  between  the  lead  and  steel  (mass,  neutron  binding 
energies,  etc.)  are  not  addressed  here. 

Flux  integrator  measurements  were  also  done  during  patient  therapy 
at  the  WPAFB  Medical  Center.  These  resulted  in  measurements  ranging 
from  1.62t0.10  mremn/radx  to  3.81±0.24  mremj^/radx  (corresponding  to  0.42 
remji  and  1.12  rem^  dose  to  the  patient),  with  a  variety  of  field  sizes 
based  on  the  particular  prescribed  treatments.  These  measurements  are 
in  the  range  of  the  measurements  done  without  patients,  and  all  measure¬ 
ments  compare  well  with  similar  studies  in  the  literature. 

Future  investigation  might  focus  on  the  calculation  of  the 
efficiency  of  the  detector (s)  used  to  count  activated  foils.  The  proce¬ 
dure  used  here  may  be  more  cumbersome  than  is  desired  or  necessary.  If 
a  gamma  counting  system  is  used  instead  of  a  beta  counting  system,  the 
dependency  on  foil  thickness  may  be  neglected  for  foils  of  similar 
thic)cness  since  gammas  are  attenuated  much  less  than  betas.  A  future 
system  might  also  be  developed  without  needing  to  use  the  standard 
graphite  pile  for  calibration:  foils  might  be  irradiated  in  the  flux 
integrators  using  a  known  neutron  flux  (as  was  done  here  in  the  reactor 
dome)  and  the  resulting  activation  compared  with  this  known  flux  to 
calibrate  the  counting  system. 

In  this  paper  the  saturation  specific  activity  5.  was  used  instead 
of  simply  the  saturation  activity  often  cited  in  the  literature.  The 
saturation  specific  activity  was  used  to  account  for  the  differences  in 
the  foils  since  their  masses  and  thic)cnesses  varied  within  a  small 
range.  Again,  if  future  work  uses  a  gamma  counter,  these  differences 
may  prove  to  be  negligible  if  foils  sufficiently  thin  are  used.  In  that 
case  there  may  be  no  need  to  account  for  the  differences,  as  was  done 
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Km  for  the  topic  of  this  thesis,  the  measurement  of  extraneous  neu¬ 
trons  from  x-ray  cancer  therapy  is  important  and  relevant  both  for  its 
own  inherent  interest  and  for  its  contribution  to  improvements  in  health 
care.  The  author  hopes  that  this  study  has,  in  some  way,  made  a  contri¬ 
bution  to  the  universal  good. 
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appendix -Jki _ Neutron  Posg  ConverilQii  Factors 


A  number  of  studiee  have  been  done  which  relate  neutron  fluence  to 
photon  dose,  neutron  dose  to  photon  dose,  and  neutron  fluence  to  neutron 
dose.  These  relations  are  usually  based  on  the  average  energy  of  the 
neutron  spectrum,  which  in  turn  depends  on  the  photon  energy  and  the 
target,  flattener,  and  collimator  assembly  materials. 

For  comparison  and  reference  a  table  of  conversion  factors  from  the 
literature  follows,  sorted  by  units.  The  references  are  given  at  the 
end  of  part  d*  For  reference  d,  the  accelerator  voltage  is  given 
instead  of  the  photon  energy  in  the  fourth  column  of  parts  a-c* 


Table  16 t  Various  Neutron  Dose  Conversion  Factors 
a.  Assorted  Units 


Factor 

Units 

Spectrum 

or 

Loc  'n 

Photon 

Energy 

Machine 

Reference 

6.03x10“5 

radn/radx 

Fast 

10  KeV 

V-C  18 

(at  400) 

5. 5x10' 

radfi/ radx 

Thermal 

10  MeV 

V-C  18 

(a:  400) 

1.2x10-5 

radn/radx 

Thermal 

18  MV 

V-C  20 

(d:  673) 

4.38 

in-beam 

16  MeV 

SL75-20 

linac 

(c:  109) 

1.5 

45  MeV 

BBC 

Betatron 

(g:  55) 

0.122 

(mrem/hr ) / 

(n/cm2-.) 

252cf 

(a:  399) 
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b.  Neutron  Fluence  per  Photon  Rad 


Factor 

Unite 

Spectrum 

or 

Loc’n 

Photon 

Energy 

Machine 

Reference 

8x10^ 

n/cm^-radx 

Faet 

18  MV 

Varian 

linac 

{d:  673) 

1.2x103 

Thermal 

10  MeV 

V-C  18 

(a:  400) 

1.8x10^ 

n/cm^-radx 

Fast 

33  MeV 

BBC 

Betatron 

(e:  400) 

1.8x105 

Past 

18  MV 

V-C  20 

(d:  673) 

2.7x10^ 

n/cxn^-radx 

Thermal 

18  MV 

V-C  20 

(d:  673) 

1.52x10^ 

n/cm^-radx 

Fast 

5  cm  out 

10  MeV 

V-C  18 

(a:  400) 

1.6xl0l0 

n/-radx 

infinite 

target 

25  MeV 

(b:  82) 

4.2x108 

n/-radx 

infinite 

target 

10  MeV 

V-C  18 

(b:  83) 

2.1x105 

n/cm2-radx 

25  NeV 

(b:  84) 

l.SxloS 

n/cm2-radx 

out  of  beam 

16  MeV 

SL75-20 

linac 

(c:  116) 

7.2x104 

n/cm^-radx 

out  of  beam 

16  MeV 

MEL  linac 

(e:  400) 

1.7x105 

n/cm2-radx 

5  cm  out 
of  beam 

19  MeV 

Siemens 

Betatron 

(e:  400) 

c.  Rad  Neutron  Dose  per  Neutron  Fluence 


Factor 

Units 

Spectrum 

Photon 

Energy 

Machine 

Reference 

3.97x10"9 

rad^-cm^/n 

252c£ 

(a:  399) 

4.57x10-10 

radn~cm2/n 

Thermal 

10  MeV 

V-C  18 

(a:  399 

4.4x10-10 

radn"cm2/n 

Thermal 

18  MV 

V-C  20 

(d:  673) 

3.0x10-9 

radn-cm^/n 

252cf 

(e:  400) 

2.74x10-11 

radn-cm^/n 

Thermal 

(f:  467) 
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d.  Hem  Neutron  Dose  per  Neutron  Flux 


Factor 

Units 

Spectrum 

Photon 

Energy 

Material 

Reference 

4.2x10“^ 

n-remn/cm^ 

Tungsten 

(b:  77) 

2.4x107 

n-remn/cm2 

E  •  2  MeV 

Lead 

(b:  77) 

V*C  refers  to  Varian  Clinac  machines • 

References:  a:  (McGinley  et  al,  1976) 

b:  (McCall  and  Swanson,  1979) 
c:  (Axton  and  Bardell,  1975) 
d:  (Purdy  and  Glasgow,  1984) 
e:  (Wilenzick  et  al,  1973) 
f:  (Attix,  1986) 
g:  (Rossano,  1989) 

From  the  last  two  entries  of  part  b.  of  the  preceding  table,  it  is 
reasonable  to  expect  between  7.2 x  10^  n/cm^-radx  and  1.7x  lO®  n/cm^-radx 
for  the  Varian  Clinac  1800*  Based  on  an  average  neutron  energy  of 
approximately  1  MeV  for  tungsten  photoneutrons  (Toms  and  Stephens,  1957: 
77-81)  a  conversion  factor  of  2.38  x  lO’®  rem-n/cm^  was  used  for  neutron 
flux  to  dose  calculations  for  the  Varian  Clinac  1800  (McCall  and  Swan¬ 
son,  1979:  77).  If  the  predominate  material  in  the  accelerator  head 
were  lead  (as  with  the  CGM  Saturne  1  accelerator  at  Kettering),  the 
average  photoneutron  energy  would  be  about  2  MeV  (NCRP  79,  1979:  30)  and 
a  factor  of  4.17xi0‘®  rem-n/cm^  would  be  appropriate  for  the  neutron 
flux  to  neutron  dose  conversion  factor.  (NCRP  79,  1979:  45;  McCall  and 
Swanson,  1979:  77).  Note  these  two  conversion  factors  are  the  inverse 
of  the  factors  in  part  d.  of  the  above  table. 
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Once  the  neutron  fluence  ie  knovm,  then,  the  neutron  dose  can  be 
calculated  by  the  formula 

// »  (2.38X  lO‘“rem,pn/cm^)<J>  (16) 

for  the  neutrone  from  the  WPAFB  electron  accelerator,  which  uses  tung¬ 
sten.  At  Kettering  Memorial  Hospital  the  formula  is 

H-(4.17x  10’®rcm„‘n/cm^)4>  (17) 

since  the  CGM  Saturne  1  uses  lead. 
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Appendix  Bi  Foil  Paramatera 


I 


The  list  which  follows  describes  each  of  the  foils  used  for  this 
report  * 


Foil 

Mass 

0i2uneter 

Thickness 

Number 

(mg) 

[cm] 

(mg/cm2 ) 

1 

125.1 

1.281 

97.1 

2 

131.1 

1.277 

102.4 

3 

134.9 

1.287 

103.4 

4 

130.1 

1.288 

99.9 

5 

135.9 

1.279 

105.8 

6 

135.3 

1.288 

103.8 

7 

132.2 

1.288 

101.5 

8 

128.8 

1.288 

98.9 

15 

136.9 

1.280 

106.4 

16 

127.1 

1.290 

97.2 

17 

125.0 

1.282 

96.8 

18 

124.8 

1.278 

97.3 

19 

128.9 

1.281 

100.0 

20 

131.9 

1.277 

103.0 

21 

133.8 

1.276 

104.6 

22 

254.0 

2.555 

49.5 

23 

242.9 

2.555 

47.4 

24 

253.2 

2.547 

49.7 

25 

246.0 

2.552 

48.1 

26 

249.2 

2.553 

48.7 

27 

255.5 

2.552 

50.0 

28 

250.6 

2.552 

49.0 

29 

240.0 

2.549 

47.0 
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Foil 

Number 

Mass 

[mg] 

30 

251.2 

31 

248.5 

32 

244.9 

33 

249.6 

34 

242.3 

35 

245.3 

36 

291.6 

37 

259.6 

38 

273.3 

39 

278.0 

40 

267.4 

41 

256.8 

42 

286.4 

43 

257.3 

.44 

287.6 

45 

265.9 

46 

278.7 

47 

271.9 

48 

271.9 

49 

300.7 

50 

279.0 

51 

281.0 

52 

266.5 

53 

266.6 

54 

266.2 

55 

266.4 

56 

276.3 

Diameter 

|cm) 

Thickness 

[mg/cm2] 

2.550 

49.2 

2.559 

48.3 

2.550 

48.0 

2.553 

48.8 

2.560 

47.1 

2.553 

47.9 

2.546 

57.3 

2.548 

50.9 

2.557 

53.2 

2.551 

54.4 

2.554 

52.2 

2.550 

50.3 

2.555 

55.9 

2.549 

50.4 

2.552 

56.2 

2.554 

51.9 

2.566 

53.9 

2.566 

52.6 

2.541 

53.6 

2.547 

59.0 

2.560 

54.2 

2.555 

54.8 

2.541 

52.6 

2.540 

52.6 

2.541 

52.5 

2.543 

52.5 

2.573 

53.1 
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Foil 

Number 

Mass 

[mgj 

57 

299*7 

58 

264.1 

59 

273.4 

60 

273.3 

61 

309.5 

62 

262.8 

63 

270.2 

64 

275.0 

65 

306.0 

66 

298.1 

67 

310.0 

68 

310.5 

69 

309.5 

70 

271.9 

71 

276.5 

72 

278.9 

73 

272.2 

74 

301.8 

75 

272.2 

76 

267.1 

77 

278.4 

78 

267.3 

79 

285.9 

80 

287.1 

61 

266.1 

Diameter 

[cm) 

Thickness 
[mg/cm2 ) 

2.575 

57.5 

2.593 

50.0 

2.571 

52.7 

2.580 

52.3 

2.578 

59.3 

2.594 

49.7 

2.562 

52.4 

2.577 

52.7 

2.583 

58.4 

2.577 

57.2 

2.584 

59.1 

2.576 

59.6 

2.586 

58.9 

2.582 

51.9 

2.579 

52.9 

2.575 

53.6 

2.559 

52.9 

2.574 

58.0 

2.571 

52.4 

2.584 

50.9 

2.581 

53.2 

2.550 

52.3 

2.555 

55.8 

2.551 

56.2 

2.547 

52.2 
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Foil 

Number 

Mass 

Img] 

Diameter 

[cm] 

Thickness 
[mg/cm^  J 

82 

298,0 

2.540 

58.8 

83 

294,5 

2.541 

58.1 

84 

263.1 

2.561 

51.1 

J-01 

119.5 

2.558 

23.3 

J-02 

122.9 

2.529 

24.5 

J-03 

122.9 

2.560 

23.8 

J-04 

481.5 

2.493 

98.6 

J-05 

496.5 

2.531 

98.7 

J-06 

468.1 

2.490 

96.1 
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i^ppendix  C;  Increase  in  MS80  and  M1170  Outputs 


The  equation  for  computing  the  output  of  a  Pu-Be  source  is 


■^a^X^Nj- - ^(e  =*  -e  M 


(29) 


where  a  is  the  average  neutron  yield  per  alpha  particle,  K  is  the  disin¬ 
tegration  constant,  N  is  the  number  of  atoms  at  time  t»0,  and  the  sub¬ 
scripts  1,  2,  3,  and  4  refer  to  239pu^  240pu^  241pu^  and  241;^ 
respectively  (Anderson,  1968:142) •  For  most  such  neutron  sources  -  as 
is  true  of  M580  and  M1170,  used  in  this  thesis  project  -  the  amount  of 
24lAm  is  negligible  or  non-existent  at  time  t*0;  thus  N4=0  and  the  third 
term  above  goes  to  zero.  Also  the  decay  of  239pu  and  240pu  ie  very 
small  so  that  the  exponentials  in  the  first  two  terms  are  approximately 
unity.  They  can  be  combined  into 

Q(0)-ai\iA^i  +  02X2^2  (30) 


which  is  the  initial  emission  rate  given  for  the  source  (Anderson, 
1968:143).  With  these  simplifications  £q  (29)  above  may  be  written  as 

A3  A4 

Now  A3*  4.83x10"’2  yr”i  and  X4  *  1.60x10“^  yr"^  (Chart  of  Nuclides, 
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1983:47)  and  a4*S2xlO”®  neutrons/a  (Mound,  1962:F-7).  Also  N3  may  be 
written  as 


where  M  ie  the  total  maoB  of  Pu  in  the  source,  P3  ia  the  fraction  of  the 
mass  which  is  241pu^  j^a  Avogadro's  number,  and  AW3  is  the  atomic 
weight  of  241pu,  Substituting  these  values  and  making  appropriate  con-- 
versions  (such  as  years  to  seconds)  gives 

Q(O  =  <?(0)-^6,84A:10‘Mp3(e‘’'’‘-e’^‘'')  (33) 


in  units  of  neutrons /second. 

For  neutron  source  M580,  the  variables  have  the  following  values: 

Q(0)  *  8.86x10^  n/s  on  23  June  1960 
t  *  29.19  years  on  1  September  1989 
M  *  79.79  g 
P3  »  0.0044 

(Mound,  1960:1;  Anderson,  1967:145;  1980:429).  Substituting  these  val¬ 
ues  into  Eq  (33)  gives  Q(t)  *  1.06x10*^  n/s  as  the  output  during  the 
time  of  this  thesis. 

(A  note  on  P3  for  M580  is  in  order  here.  The  shipping  data  from 
Monsanto  for  M580  could  not  be  found,  so  a  specific  reference  for  p3  was 
not  available  to  this  researcher.  Presumably  Rossano  had  the  reference 
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but  it  has  since  been  lost.  There  are  at  least  two  ways  to  determine 
the  value  for  p3  for  the  time  of  this  study,  then.  Using  the  value 
reported  by  Rossano  for  MS80  at  t-26.7  years,  one  could  solve  Eg  (33) 
^or  P3  since  every  other  value  was  known.  At  t=26.7  years,  Rossano 
reported  Q(t )=10. 5x10^  n/s  (Rossano,  1969:64).  Solving  Eg  (33)  for  P3 
gives  0.0044  .  Another  method  is  to  compare  the  current  output  of  M560 
and  M1I70  to  each  other,  as  is  done  in  Appendix  G.  The  ratios  of 
measured  output  compared  with  the  ratios  of  calculated  output  from  Eg 
(33)  allows  one  to  solve  for  P3  for  M580,  with  P3  well  known  for  M1170 
but  not  BO  well  known  for  M560.  In  addition  to  these  calculations  there 
are  two  literature  references  -  (Anderson,  1967:145;  1980:429)  -  which 
lend  support  to  this  induced  value  of  P3  for  MS80. ) 

Since  M580  is  the  neutron  source  for  the  standard  graphite  pile,  an 
increase  in  its  flux  causes  an  increase  in  the  thermal  neutron  flux  of 
the  pile.  Dividing  Q(t)  by  Q(0)  gives  the  ratio  Q(t)/Q(0)  =  1.193  , 

hence  the  assertion  that  the  neutron  output  of  PuBe  source  H560  has 
increased  19.3%  since  its  birth  (where  t=29.19  years). 

For  neutron  source  M1170,  the  following  values  were  used: 

Q(0)  =  9.00x106  n/B  on  2  March  1962 
t  =  27.38  years  on  20  July  1989 
M  -  76.36  g 
P3  =  0.0055 

(Mound,  1962 :F-1).  These  values  in  Eg  (33)  give  Q(t)  -  1.1x10*^  n/s 
for  the  output  during  the  time  of  this  thesis* 
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Appandir  D!  Foil  Thickne»«  v«  Count »  Data 


The  two  table  in  this  appendix  give  the  data  used  for  foil  thick¬ 
ness  dependency  determination.  Table  17  gives  the  data  used  in  Figure 
3,  for  saturation  counts/thickness  vs  foil  thickness.  Table  16  gives 
the  data  used  in  Figure  4,  the  graph  of  saturation  counts/mass  vs  foil 
thickness.  The  saturation  counts  (units  were  computed  by 


(1  -g''"'') 


(18) 


Thickness  is  in  units  of  mg/cm^  ,  and  mass  is  in  units  of  milligram. 


Table  17;  Data  for  Foil  Thickness  vs  Saturation  Counts/Thickness 


Foil 

Count  No. 

Thickness, X 
(mg/cm2 ) 

Cfl/x 

(cm^/mg-B) 

253 

24.8 

2.015  ±  0.020 

254 

24.8 

1.995  ±  0.020 

255 

24.8 

1.987  ±  0.020 

256 

24.6 

2.007  ±  0.020 

257 

25 

1.991  ±  0.020 

258 

25 

2.021  ±  0.022 

259 

25 

2.021  ±  0.020 

260 

24.9 

2.002  1  0.020 

261 

24.9 

2.023  ±  0.024 

262 

24.9 

1.977  ±  0.021 

263 

24.6 

2.058  ±  0.021 

264 

24.6 

2.049  ±  0.021 

265 

24.6 

2.048  ±  0.021 

266 

49.6 

1.515  ±  0.013 

267 

49.6 

1.536  ±  0.014 

268 

74.7 

1.1857  ±  0.0092 

269 

74.7 

1.217  ±  0.011 

272 

49.5 

1.518  ±  0.013 

273 

49.5 

1.541  ±  0.015 

274 

74.5 

1.1949  ±  0.0093 

275 

74.5 

1.218  ±  0.011 
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Table  18:  Data  for  Foil  Thickness  vs  Saturated  Counts/Mass 


Foil 

Count  No. 

Thickness, X 
[mg/cm2 ] 

Cg/X 

(cm^/mg-s) 

253 

24.8 

1.542  ±  0.015 

254 

24.8 

1.527  ±  0.015 

255 

24.8 

1.521  ±  0.015 

256 

24.8 

1.536  ±  0.015 

257 

25 

1.518  ±  0.015 

258 

25 

1.S40  ±  0.016 

2S9 

25 

1.540  ±  0.016 

260 

24.9 

1.534  ±  0.016 

261 

24.9 

1.550  ±  0.018 

262 

24.9 

1.515  ±  0.016 

263 

24.6 

1.549  t  0.016 

264 

24.6 

1.541  ±  0.016 

265 

24.6 

1.541  ±  0.016 

266 

49.8 

1.1572  t  0.0097 

267 

49.8 

1.175  ±  0.011 

268 

74.7 

0.9066  ±  0.0071 

269 

74.7 

0.9303  ±  0.0081 

272 

49.5 

1.1523  i  0.0097 

273 

49.5 

1.170  ±  0.011 

274 

74.5 

0.9084  ±  0.0070 

275 

74.5 

0.9261  ±  0.0081 

*^’0  show  how  these  numbers  were  obtained,  take  Foil  Count  No.  253  in 
tabic  17.  The  values  for  each  of  the  variables  in  Eq  (18)  are  taken 
from  Appendix  I,  which  contains  the  data  from  each  irradiation  and 
counv.*  For  #253  the  following  values  were  used: 

tv/  s  8  minutes 
te  »  2875  minutes 
tc  »  500  seconds 

Cnet  *  10175  -  (500  s) (0.320  background)  *  10015  counts 
ti,  «  2.6935  days 

X  *  (ln2/2.6935  days)  *  (1  day/86400  seconds) 
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These  values  were  substituted  into  the  equation;  the  value  of  Cg 
thus  obtained  was  divided  by  the  value  of  the  thickness  in  the  second 
column,  to  give  the  value  in  the  third  column.  All  following  values, 
and  the  values  in  the  next  table,  were  calculated  similarly. 
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^pendi%  Et  Detector  Efficieqey 


The  data  for  the  graph  in  Figure  6  follow.  This  data  was  used  to 
determine  the  constant  fx  in  Eq  (23)  for  detector  efficiency,  €. 


Table  19:  Cadmium  Difference  Pairs  in  Graphite  Pile 


Ser . 

No. 

Str 

B/C 

C,  *  0 

e 

*±a 

♦  /« 

<t> 

16 

B 

b 

110.1  t  1.2 

0.5814 

± 

0.0018 

1658  ± 

33 

3854 

19 

B 

c 

41.00  t  0.51 

0.5694 

± 

0.0018 

64 

B 

b 

107.7  ±  1.4 

0.5961 

± 

0.0018 

1673  ± 

40 

3854 

82 

B 

c 

38.04  ±  0.68 

0.6037 

± 

0.0019 

13 

2 

b 

235.8  ±  1.7 

0.5558 

± 

0.0017 

3239  ± 

50 

8256 

22 

2 

c 

91.82  ±  0.97 

0.5771 

± 

0.0018 

67 

2 

b 

224.6  ±  2.0 

0.6094 

± 

0.0019 

3489  ± 

60 

8256 

85 

2 

c 

83.3  ±  1.0 

0.6056 

± 

0.0019 

10 

3 

b 

249.0  ±  1.7 

0.S9C6 

i 

0.0018 

4050  ± 

53 

9449 

25 

3 

c 

82.9  ±  1.0 

0.5711 

± 

0.0018 

70 

3 

b 

245.4  ±  2.1 

0.5949 

± 

0.0018 

4084  i 

61 

9449 

88 

3 

c 

77.38  ±  0.97 

0.6005 

± 

0.0019 

7 

n 

b 

155.7  ±  1.3 

0.5576 

± 

0.0017 

2870  ± 

34 

6824 

28 

B 

c 

32.37  ±  0.45 

0.5398 

± 

0.0017 

73 

n 

b 

153.2  ±  1.6 

0.6049 

± 

0.0019 

3125  ± 

45 

6824 

91 

n 

c 

27.59  ±  0.58 

0.6113 

± 

0.0019 

n 

5 

b 

72.83  ±  0.93 

0.5912 

t 

0.0018 

1647  ± 

24 

3794 

i 

El 

5 

c 

5.98  ±  0.18 

0.5676 

± 

0.0018 

1 

76 

5 

b 

74.29  t  0.95 

0.6011 

± 

0.0019 

1696  ± 

25 

3794 

1 

94 

5 

c 

6.04  i  0.12 

0.6062 

± 

0.0019 

1 

6 

D 

33.85  ±  0.63 

0.5796 

± 

0.0018 

793  i 

16 

1754 

34 

6 

B 

0.96  ±  0.10 

0.5640 

± 

0.0017 

! 

79 

6 

b 

34.48  ±  0.65 

0.5930 

± 

0.0018 

825  ± 

16 

1754 

97 

6 

c 

1.052  ±  0.041 

0.5495 

± 

0.0017 

For  the  above  table,  Cg  was  calculated  from  Eq  (18);  €*  was  calcu¬ 
lated  from  the  thickness  x  with  the  expression  given  on  page  29; 
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was  calculated  by  re-arrangement  of  £q  (23);  and  was  taken  as  1*193 
times  the  original  thermal  flux  in  the  standard  graphite  pile  (see 
Appendix  C) •  For  clarity,  the  pertinent  equations  are  reproduced  below: 


C. 


X.  O  6 


(1  -c*’''‘)(l 


(18) 


e’  =  exp[(-0.01045)x] 


<!>,*/* 


1.128/11/^^^^  C,ib)  FcgC.CCd) 

0  0^'  a  _  b  Cd 


(23a) 


A  linear  least  squares  of  ^thfx  vs  from  the  above  table  gave 
f^a(0.428±0.021)  ,  to  complete  the  efficiency  term  of  the  detector 

system* 
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The  data  used  for  Table  2,  Neutron  Fluxes  in  Graphite  Pile,  are 
given  here.  Equation  (23)  was  used  to  calculate  the  fluxes,  using  the 


value  for  fx  determined  from  Appendix  E.  The  values  used  for  these 
calculations  are,  as  for  all  pertinent  calculations,  listed  by  count 
number  in  Appendix  I.  The  values  for  in  Table  2  are  averages  of  the 
values  below  for  each  stringer  location. 


appendix  fl I - CQaoarigQn  of  MSflO  and  M117Q  Outputs 


This  appendix  describes  the  measures  used  to  compare  the  output  of 
neutron  sources  Ml 170  and  M580.  Their  outputs  were  both  measured  on  the 
same  day.  The  setup  was  as  in  Figure  9.  A  helium-3  tube  with  a  2-inch 
cadmium  cap  to  block  direct  neutrons  was  used  in  the  long  counter.  The 
long  counter  was  set  on  the  edge  of  the  Ping-pong  table  in  the  basement 
of  Bldg  470.  Parameters  of  the  experiment  are  listed  in  Table  20. 


Figure  9  Setup  for  Comparison  of  M580  and  Ml 170  Neutron  Sources 


Four  counts  of  100  seconds  each  were  made  with  both  sources*  The 
results  are  listed  in  Table  21* 
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Table  20:  MS80  and  M1170  Neutron  Source  Comparison  Parameters: 


Distance,  source  to  He-3  tube: 

2  meters 

Height  of  source  and  counter: 

94  cm 

Long  counter  settings: 

1200  volts 

Table  21:  M580  and  Ml 170  Comparison 


Trial 

M580  Counts 

Ml  170  Counts 

1 

34271 

•  36459 

2 

33932 

36422 

3 

34454 

36471 

4 

34244 

36419 

Average : 

34225  ±  217 

36443  ±  26 

The  ratio  of  measured  outputs  is  1.0648  ±  0.0068  .  The  measure¬ 
ments  were  done  on  31  Oct  89.  On  this  date,  t  »  26.67  yrs  for  M1170 
and  t  *  29.36  yrs  for  M580.  The  calculated  values  according  to  Eq. 
(33)  are  1.1x10*^  n/s  for  M1170  and  1.06x10*^  n/s  for  M580  (using 
P3»0.0044  ).  The  calculated  ratio,  then,  is  1*04*  This  agrees  reason¬ 

ably  well  with  the  measured  ratio,  and  confirms  0*0044  as  the  original 
241py  fraction  in  M580* 
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Appendix  Hs  Flux  Integrator  Calibration  Congtant 


The  value  of  the  calibration  constant  K  for  the  flux  integrators^ 
used  in  the  equation  ,  was  obtained  by  linear  least  squares. 

The  graph  of  the  data  was  presented  in  section  4.3.  This  appendix 
describes  the  data  and  calculations  used  to  obtain  the  calibration  con¬ 
stant. 

To  calculate  this  value,  one  may  let  4>  be  the  independent  variable, 
or  one  may  let  5.  be  the  independent  variable.  Here  5.  was  plotted  on 
the  x-axis  and  4>  along  the  y-axis  to  match  the  equation  above.  This 
assigns  5.  as  the  independent  value. 

5.  was  calculated  as  /I. /Mass,  with  A,  from  Eq.  (7).  4)  was  deter¬ 

mined  by 


(34) 


where  Q(M1170)  is  the  output  of  PuBe  source  M1170,  1.15  is  the  aniso¬ 
tropy  of  M1170  (Rossano,  1969:  69),  and  ground  scatter  was  accounted  for 
by  Eq  (26) . 

The  method  of  linear  least  squares  was  used.  The  equations  after 
the  table  are  from  chapter  six  of  Bevington.  The  following  table  gives 
the  data  and  the  statistical  expressions  used. 
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Table  22:  Flux  Integrator  Calibration,  I 


Cnt 

# 

■ 

iS.±a) 

X  10-6 

[Bq/mg] 

n 

<)>' 

4)  S.  ^  0 

a  6  X', 

n 

302 

B 

3.44 
±  0.14 

104 

10816 

3.58 
±  0.15 

1.183 
±  0.096 

109.39 

303 

2 

4.91 
±  0.19 

151 

22801 

7.41 
±  0.29 

2.41 
±  0.10 

153.04 

4.16 

304 

3 

4.80 
±  0.19 

151 

22801 

7.25 
±  0.29 

2.30 
±  0.18 

149.77 

Bi 

305 

B 

3.23 

±0.13 

104 

10816 

3.36 

±0.14 

1.043 
±  0.084 

103.16 

0.71 

306 

5 

3.29 

±0.13 

104 

10816 

3.42 
±  0.14 

1.082 
±  0.085 

104.94 

0.88 

332 

6 

3.82 
±  0.15 

121 

14641 

4.62 
±  0.18 

1.46 
±  0.11 

120.68 

0.10 

333 

B 

3.72 
±  0.15 

121 

14641 

4.50 
±  0.18 

1.38 
±  0.11 

117.71 

10.82 

334 

8 

3.42 
±  0.14 

107 

11449 

3.66 
±  0.15 

1.170 
±  0.096 

108.80 

3.24 

335 

9 

3.21 

±  0.13 

107 

11449 

3.43 
±  0.14 

1.030 
±  0.042 

102.56 

19.71 

Sum 

N=9 

33.84 
±  0.46 

1070 

130230 

41.23 
±  0.58 

13.06 
±  0.36 

70.18 

a  -  — - — - — - -  ■  ■  -  7.26 


NUs.-ys.y4> 

5 - i= - _  2968.82 


‘0  026 


0!- 


al  = 


A 
Ns^ 


-  43.28 


29825 


so  that 


so  that 


a„  -  6.58 


0,-  173 
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These  equations  give 

=  (2970±  170)S,  +  (7.3±6.6)  (28) 

for  the  relation  between  flux  and  saturation  activity.  This  is  in 
agreement  with  Rossano's  results,  though  the  uncertainties  are  larger 
than  his. 

The  uncertainties  of  the  data  are  given  in  the  table,  but  they  are 
not  used  to  calculate  the  final  uncertainly  in  the  activity-to-f lux 
equation.  The  uncertainties  were  first  calulated  using  propogation  of 
error,  as  outlined  in  Bevington.  However,  the  uncertainties  obtained 
were  not  believable  as  calculated;  the  uncertainty  in  K  was  found  to  be 
about  3/4  of  the  value  of  K  itself.  If  one  were  to  calculate  K  for  each 
data  point,  the  individual  uncertainties  were  never  more  than  about  5% 
of  the  value  of  K.  Either  the  uncertainties  by  propogation  of  errors 
blew  up  by  this  stage  in  the  calculations,  or  they  were  misapplied.  At 
any  rate  it  was  decided  to  use  the  uncertainties  given  here,  based  on 
the  statistical  method  used  to  calculate  K  from  4>  and  S..  They  gave 
much  more  believable  results. 

If  4>  were  used  as  the  independent  variable,  so  that  4>  is  plotted 
along  the  x-axis  and  along  the  y-axis,  the  last  two  columns  of  the 

table  above  and  all  the  equations  that  follow  will  change.  For  compari¬ 
son,  they  are  given  here. 
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Table  23:  Flux  Integrator  Calibration,  II 


Cnt  # 

i 

a  •»'  6  Xj 

Af 

302 

1 

0.0327 

2.89x10-6 

303 

2 

0.0482 

8.1x10-7 

304 

3 

0.0482 

4.0x10-8 

305 

4 

0.0327 

1.6x10-7 

306 

5 

0.0327 

4.0x10-8 

332 

6 

0.0383 

1.0x10-8 

333 

7 

0.0383 

1.21x10-6 

334 

8 

0.0337 

2.5x10-7 

335 

9 

0.0337 

2.56x10-6 

Done  this  way,  the  etjuations  give  these  results: 


A  =  27170 


a  - 


Y.siz^-Y_s.Y.^s. 


--0.0017 


nY.^S.-Y,S.Y.^ 

b - ^ - ^-3.30585x10'^ 


1*39x10- 


^Isl 


al - 5.46x10'* 

oS  —  =  3.77x10'^ 


SO  that 


so  that 


a„  -  2.34>:10' 


0^*  1.94x10’ 


Under  this  arrangement,  A' -  1 /5  «  3025*  178  ,  if  the  appropriate 
terms  are  inverted.  The  uncertainty  in  the  intercept  is  greater  than 
the  intercept  itself,  so  again  the  intercept  may  be  interpreted  as  zero. 
Statistically  the  answer  is  not  different  from  the  previous  one;  and 
based  on  the  experiment  used  <(>  may  be  more  the  independent  variable  than 
5.  .  This  arrangement  was  used  by  a  previous  researcher  (Rossano,  1988: 
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37-39).  However,  this  researcher  chose  to  use  the  first  arrangement 
given  above  in  his  calculations;  it  seemed  to  be  the  more  direct  method 
to  reach  the  desired  equation.  Incidentally  it  also  gives  the  same 
answer  as  in  the  previous  work  on  this  topic  (Rossano,  1988:  38),  though 
it  was  not  chosen  specifically  for  this  reason. 
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Appendix  I:  Foil  Counting  Data 


The  following  list  includes  most  of  the  pertinent  irradiation  data 
from  each  foil  measurement  done  in  preparation  for  this  thesis.  The 
count  number  is  simply  an  arbitrary  designator  for  each  individual  count 
of  irradiated  foils.  The  foil  numbers  are  those  given  in  the  previous 
section.  C(tc)  refers  to  the  gross  (or  total)  counts  in  counting  time 
t^;  no  correction  for  background  counts  are  included  in  that  column,  tg 
is  the  foil  exposure  time,  and  t^^  is  the  waiting  time  between  end  of 
irradiation  and  beginning  of  counting,  r^g  is  the  background  count  rate 
for  the  particular  detector  used,  as  measured  by  background  count  data 


for  some  time 

within 

i  a  few  days 

of  the  particular 

count . 

All  counting 

was  done  with 

the  G- 

M  detector  with  a 

silver  backscatterer 

It  was 

located  in  Building 

470,  in  the 

lab  rooms. 

Count 

Foil 

C(tc) 

^c 

^bg 

te 

^w 

Number 

Number 

iB] 

(s-Jj 

(m] 

[m] 

1 

40 

3110 

400 

0.2772 

1404 

7 

4 

41 

6287 

400 

0.2772 

1345 

40 

7 

42 

13815 

400 

0.2772 

1396 

15 

10 

43 

20692 

400 

0.2772 

1301 

19 

13 

44 

20082 

400 

0.2772 

1340 

29 

16 

45 

9244 

400 

0.2772 

1302 

3 

19 

46 

6987 

800 

0.2772 

1308 

52 

22 

47 

9323 

500 

0.2772 

1255 

20 

25 

48 

6751 

400 

0.2772 

1254 

9 

28 

49 

5578 

800 

0.2772 

1307 

32 

31 

50 

1582 

1000 

0.2772 

1380 

4 
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Count 

Number 

Foil 

Number 

C{tc) 

(S) 

fbg 

[8-2] 

te 

im] 

[»] 

34 

51 

485 

1000 

0.2772 

1381 

23 

36 

37 

7535 

500 

0.2772 

1370 

20 

64 

22 

6371 

400 

0.2772 

2929 

5765 

67 

23 

13170 

400 

0.2772 

2933 

5781 

70 

24 

14512 

400 

0.2772 

2968 

5749 

73 

25 

9145 

400 

0.2772 

2967 

5743 

76 

26 

6487 

600 

0.2772 

2895 

5852 

79 

27 

3096 

600 

0.2772 

2895 

5865 

82 

31 

3164 

100 

0.2772 

9917 

37 

85 

32 

6916 

100 

0.2772 

9919 

24 

88 

33 

6397 

100 

0.2772 

9920 

47 

91 

34 

2303 

100 

0.2772 

9922 

37 

94 

35 

2622 

500 

0.2772 

9923 

48 

97 

36 

1139 

1000 

0.2772 

9923 

68 

141 

36 

2566 

1000 

0.284 

0.83 

214  ± 

10 

142 

38 

1974 

1000 

0.284 

0.83 

233  t 

10 

143 

42 

2048 

1000 

0.284 

0.85 

245  ± 

10 

144 

40 

1934 

1000 

0.284 

0.85 

285  ± 

10 

145 

35 

2126 

1000 

0.284 

0.83 

273  t 

10 

146 

45 

2092 

1000 

0.284 

0.83 

283  ± 

10 

147 

39 

1500 

1000 

0.284 

0.83 

311  ± 

10 

148 

44 

2038 

1000 

0.284 

0.83 

320  ± 

10 

149 

46 

3228 

1000 

0.303 

0.92 

505 

150 

47 

1412 

1000 

0.303 

0.64 

503 

151 

48 

2152 

1000 

0.303 

2.17 

447 
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Count 

Foil 

C(tc) 

tc 

[8-1) 

te 

Number 

Number 

[S] 

[m] 

[m] 

152 

49 

3151 

1000 

0.303 

0.80 

443 

253 

J-08 

10175 

500 

0.320 

2875 

8 

254 

J-08 

10056 

500 

0.320 

2875 

20 

255 

J-08 

9996 

500 

0.320 

2875 

29 

256 

J-08 

10082 

500 

0.320 

2875 

40 

257 

J-09 

10062 

500 

0.320 

2875 

51 

258 

J-09 

9174 

600 

0.320 

2875 

1699 

259 

J-09 

10192 

500 

0.320 

2875 

61 

260 

J-10 

10041 

500 

0.320 

2875 

72 

261 

J-10 

7652 

500 

0.320 

2875 

1678 

262 

J-10 

8961 

600 

0.320 

2875 

1688 

263 

J-11 

9874 

500 

0.320 

2875 

255 

264 

J-11 

9811 

500 

0.320 

2875 

264 

265 

J-11 

9789 

500 

0.320 

2875 

274 

266 

J-08, 

J-09 

14553 

500 

0.320 

2875 

285 

267 

J-08, 

J-09 

11481 

500 

0.320 

2875 

1712 

266 

J-10, 

J-09, 

J-08 

17021 

500 

0.320 

2875 

297 

269 

O  00 

o  o 

1  1  1 

13567 

500 

0.320 

2875 

1724 

270 

J-11, 

18973 

500 

0.320 

2875 

311 

J-10, 

J-09, 

J-08 


80 


Count 

Foil 

C(tc) 

tc 

rb« 

[B-i] 

Number 

Number 

to] 

[m] 

tm) 

271 

J-11, 

J-10, 

J-09, 

J-08 

19017 

500 

0.320 

2875 

320 

272 

J-10, 

J-11 

14375 

500 

0.320 

2875 

331 

273 

J-10, 

J-11 

11352 

500 

0.320 

2875 

1755 

274 

J-09, 

J-10, 

J-11 

16974 

500 

0.320 

2875 

341 

275 

J-09, 

J-10, 

J-11 

13500 

500 

0.320 

2875 

1744 

276 

J-08, 

J-09, 

J-10, 

J-11 

18895 

500 

0.320 

2875 

351 

277 

J-08, 

J-09, 

J-10, 

J-11 

14930 

500 

0.320 

2875 

1735 

298 

26 

8214 

4000 

0.319 

1.78 

165 

299 

31 

10537 

5000 

0.319 

0.75 

163 

300 

33 

6985 

5000 

0.319 

0.82 

243 

301 

34 

8754 

6000 

0.319 

0.68 

75  t  5 

302 

25 

10654 

10000 

0.321 

2382 

29 

303 

37 

11831 

10000 

0.3219 

1804 

56 

304 

43 

11359 

10000 

0.3219 

1804 

226 

305 

35 

13370 

10000 

0.325 

3891 

18 

306 

36 

13997 

10000 

0.325 

3891 

208 

307 

38 

12037 

6000 

0.324 

0.63 

107 

308 

39 

12298 

8000 

0.324 

0.83 

214 
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Count 

Foil 

C(tc) 

Number 

Number 

[®1 

[m] 

[m] 

309 

40 

12226 

8000 

0.324 

0.70 

332  t  5 

310 

51 

2644 

1000 

0.335 

1 

313 

311 

32 

10774 

5000 

0.335 

1 

1477 

312 

23 

11326 

5000 

0.335 

1 

1388 

313 

24 

2779 

1000 

0.335 

1 

323 

314 

27 

12961 

5000 

0.335 

1 

660 

315 

30 

9246 

4000 

0.335 

1 

1285 

316 

44 

2671 

1000 

0.335 

1 

326 

317 

42 

9919 

4000 

0.335 

1 

1034 

318 

53 

2858 

1000 

0.335 

1 

334 

319 

45 

10185 

4000 

0.335 

1 

1079 

320 

55 

14959 

1000 

0.335 

6 

312 

321 

54 

6350 

1000 

0.335 

6 

343 

322 

52 

5251 

1000 

0.335 

6 

363 

323 

56 

3032 

1000 

0.335 

6 

382 

324 

58 

14624 

1000 

0.335 

6 

470 

325 

57 

6809 

1000 

0.335 

6 

488 

326 

62 

5313 

1000 

0.335 

6 

527 

327 

63 

2423 

1000 

0.335 

6 

546 

328 

59 

15219 

1000 

0.335 

6 

373 

329 

60 

6942 

1000 

0.335 

6 

399 

330 

61 

5761 

1000 

0.335 

6 

419 

331 

64 

2589 

1000 

0.335 

6 

437 

82 


Count 

Number 

Foil 

Number 

C(tc) 

[8] 

^bg 

(B-J) 

te 

[m) 

(m] 

332 

71 

12065 

10000 

0.331 

2350 

25 

333 

74 

11439 

10000 

0.331 

2350 

474 

334 

65 

11176 

12000 

0.338 

1619 

64 

335 

66 

15533 

18000 

0.338 

1619 

269 

336 

68 

11572 

1000 

0.342 

3 

177 

337 

67 

8558 

1000 

0.342 

3 

305 

338 

69 

11731 

1200 

0.342 

3 

325 

339 

70 

7165 

1200 

0.342 

3 

248 

340 

83 

12368 

1200 

0.342 

3 

396 

341 

84 

9462 

1200 

0.342 

3 

419 

342 

72 

13481 

1200 

0.342 

3 

430 

343 

73 

14311 

1800 

0.342 

3 

452 

344 

75 

11766 

1200 

0.342 

3 

475 

345 

76 

10222 

1800 

0.342 

3 

499 

346 

79 

13349 

1200 

0.342 

3 

488 

347 

80 

14249 

1800 

0.342 

3 

510 

348 

78 

11037 

1200 

0.342 

3 

562 

349 

77 

10086 

1800 

0.342 

3 

585 

350 

46 

5946 

1200 

0.327 

6.18 

435 

351 

47 

4841 

1200 

0.327 

6.18 

457 

352 

48 

3524 

1200 

0.327 

6.18 

470 

353 

49 

5872 

1200 

0.327 

6.18 

491 

354 

40 

3792 

1200 

0.327 

6.18 

516 

355 

38 

6222 

1200 

0.327 

6.48 

538 
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Count 

Number 

Foil 

Number 

C(tc) 

CD  ft 

^bg 

(8-2) 

(m) 

[m] 

356 

39 

2682 

1800 

0,327 

6.48 

1052 

357 

41 

1432 

1800 

0.327 

6.48 

1083 

358 

25 

6215 

1800 

0.327 

6.18 

1120 

359 

22 

7219 

1800 

0.327 

6.18 

1152 

360 

50 

2507 

1800 

0.327 

6.18 

1184 

361 

24 

1477 

1800 

0.327 

6.18 

1216 

362 

26 

5945 

1200 

0.327 

6.18 

1248 

363 

31 

4179 

1200 

0.327 

6.18 

1272 

364 

51 

2378 

1800 

0.327 

6.18 

1295 

365 

44 

1365 

1800 

0.327 

6.18 

1329 

366 

33 

5246 

1200 

0.327 

6.18 

1350 

367 

34 

2868 

1200 

0.327 

6.18 

1371 

368 

27 

2307 

1800 

0.327 

6.18 

1394 

369 

42 

1409 

1800 

0.327 

6.18 

1438 

370 

37 

5776 

1200 

0.327 

6.18 

1459 

371 

• 

43 

3865 

1200 

0.327 

6.18 

1483 

372 

45 

2290 

1800 

0.327 

6.18 

1504 

373 

23 

1377 

1800 

0.327 

6.18 

1535 

374 

32 

4735 

1200 

0.327 

6.18 

1557 

375 

35 

4043 

1800 

0.327 

6.18 

1583 

376 

36 

2257 

1800 

0.327 

6.18 

1616 

377 

30 

1317 

1800 

0.327 

6.18 

1648 
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